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Abstract 

The  paper  deals  with  the  formulation  of  a  consistent  set  of  governing  equations  apt  to  de¬ 
scribe  the  physical  phenomenology  comprising  the  hypersonic  flow  field  of  an  ionized  gas 
mixture  and  the  electromagnetic  field.  The  governing  equations  of  the  flow  field  and  those 
of  the  electromagnetic  field  are  revisited  in  sequence  and  differences  or  similarities  with 
past  treatments  are  pointed  out  and  discussed.  The  equations  governing  the  flow  field  hinge 
on  the  customary  balance  of  masses,  momenta  and  energies.  The  equations  governing  the 
electromagnetic  field  are  introduced  both  directly  in  terms  of  the  Maxwell  equations  and 
by  recourse  to  the  scalar  and  vector  potentials.  The  theory  of  linear  irreversible  thermody¬ 
namics  based  on  the  entropy-balance  equation  is  also  revisited  for  the  purpose  of  obtaining, 
consistently  with  the  presence  of  the  electromagnetic  field,  the  phenomenological  relations 
required  to  bring  the  governing  equations  into  a  mathematically  closed  form.  Old  problems, 
such  as  the  influence  of  the  medium  compressibility  on  chemical-relaxation  rates  or  the  im¬ 
portance  of  cross  effects  among  generalized  fluxes  and  forces,  are  re-discussed;  additional 
problems,  such  as  the  necessity  to  consider  the  tensorial  nature  of  the  transport  properties 
because  of  the  presence  of  the  magnetic  field,  are  pointed  out.  A  non-conventional  choice 
of  first-tensorial-order  generalized  forces  and  corresponding  fluxes  is  proposed  which  ap¬ 
pears  to  offer  more  simplicity  and  better  convenience  from  a  conceptual  point  of  view  when 
compared  to  alternative  definitions  customarily  used  in  the  literature.  The  applicability  do¬ 
main  of  the  present  formulation  is  clearly  outlined  and  recommendations  for  further  work 
are  given. 
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Nomenclature 


A  vector  potential 

Ak  affinity  of  k- th  chemical  reaction 

B  magnetic  induction 

C  constant  (in  Arrhenius  law) 

c  speed  of  light  in  vacuum,  299792458  m-s_1 

D ij  diffusion  tensor 

D  ‘  thermodiffusion  tensor 

Djfc  diffusion  tensor  (Fick  law) 

dj  diffusion  vector  (kinetic  theory) 

E  electric-field  intensity 

Ea  activation  energy  (in  Arrhenius  law) 

e  electronic  charge,  1.602176462  •  10~19  C 

em  matter  energy  per  unit  total  mass 

em>v  matter-energy  production 

ecm.v  electromagnetic-energy  production 

Fi  external  force  (kinetic  theory) 

5F i  generalized  force  [Eq.  (127)] 

fi5  Helmholtz  potential  of  5-th  molecular  degree  of 
freedom  of  z-th  component 
G  generic  extensive  variable 

g  generic-variable  density  (mass) 

gv  generic-variable  density  (volume) 

g  generic- variable  production  (mass) 

gv  generic-variable  production  (volume) 

hi  enthalpy  of  i- th  component  per  its  unit  mass 

JG  generic-variable  diffusive  flux 

JEra  matter-energy  diffusive  flux 

Jm.  component-mass  diffusive  flux 

Jm*  element-mass  diffusive  flux 

j 

JQ  electric-charge  diffusive  flux  or  conduction-current  density 

Jq  heat  flux  (see  text) 

J,  internal-energy  diffusive  flux 

JUi5  diffusive  flux  of  (see  below) 

Js  entropy  diffusive  flux 

j  electric-current  density 

Kb  Boltzmann  constant,  1.3806503-10~23  J-K"1 

K‘k  chemical-equilibrium  constant  (concentrations) 

ii  number  of  molecular  degrees  of  freedom  of  /-th  component 

M  magnetization 

M  gas-mixture  average  molar  mass 

Mi  component  molar  mass 

M?  element  molar  mass 
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V 

Lri  hrb 

r^k 

Xe 

xc 


number  of  components 

Avogadro  number,  6.02214199  •  1023 

component  particle  number 

element  particle  number 

polarization 

pressure 

partial  pressure  of  ?'-th  component 

component  molar  electric  charge 

electric  charge  per  unit  mass 

universal  gas  constant,  8.314472  J-K_1 

number  of  chemical  reactions 

number  of  elements 

entropy  per  unit  total  mass 

entropy  of  i- th  component  per  its  unit  mass 

entropy  production 

entropy  production  related  to  tensorial  order  0,1,2 

temperature  (thermal  equilibrium) 

temperature  associated  with  <5-th  molecular  degree 

of  freedom  of  i- th  component 

time 

unit  tensor 

internal  energy  of  the  gas  mixture 

internal  energy  distributed  over  <5-th  molecular  degree 

of  freedom  of  i-th  component 

internal  energy  per  unit  total  mass 

internal  energy  distributed  over  5-th  molecular  degree 

of  freedom  of  i-th  component  per  unit  mass  of  (-th  component 

internal-energy  production 

production  of  Uts 

velocity  vector 

traceless  symmetric  part  of  velocity  gradient 
antisymmetric  part  of  velocity  gradient 
specific  volume 

specific  volume  of  i-th  component 
volume  production 
component  diffusion  velocity 
molar  fraction  of  (-th  component 
component  mass  fraction 
element  mass  fraction 

dielectric  constant  of  vacuum,  8.854187817  •  10-12  F-m_1 
temperature  exponent  (in  Arrhenius  law) 
reaction  constant  (forward,  backward) 
scalar  electrical  conductivity 
electrical-conductivity  tensor 
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X’-  pressoelectrical-conductivity  tensor 

X"  thermoelectrical-conductivity  tensor 

Xr  thermal-conductivity  tensor  (see  text) 

Hi  chemical  potential  of  i- th  component 

dynamic-viscosity  tensor 
uw  bulk-viscosity  coefficient 

Uki  global  stoichiometric  coefficient 

i&  v'kl  stoichiometric  coefficient  (reactant,  product) 
chemical-reaction  rate 
n  normal  mean  stress 

p  total-mass  density 

[)c  electric-charge  density 

Pi  component  partial  density 

p“  element  partial  density 

CTij  formation-matrix  coefficient 

T  stress  tensor 

Tos  traceless  symmetric  part  of  stress  tensor 

Tm  Maxwell  stress  tensor 

<X»(,  generic-variable  flux 

0  scalar  potential 

%i  generalized  force  [Eq.  (Ill)] 
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1  Introduction 


Interest  in  theoretical  investigations  [1-26]  aimed  at  the  understanding  of  the  fluid 
dynamics  and  the  thermodynamics  of  flows  subjected  to  the  action  of  electric  and/or 
magnetic  fields  can  be  traced  back,  at  least,  to  the  first  half  of  the  past  century.  The 
main  driving  motivation  that  justifies  such  an  interest  was  probably  best  expressed 
by  Resler  and  Sears  [9]  in  1958: 

If  a  fluid  is  a  conductor  of  electricity,  the  possibility  arises  that  an  electric  body 
force  may  be  produced  in  it  that  will  affect  the  fluid  flow  pattern  in  a  significant 
way  . . .  The  attractive  thing  about  the  electric  body  force  ...  is  that  it  can  be 
controlled,  insofar  as  the  current  and  the  magnetic  field  can  be  controlled,  and 
perhaps  made  to  serve  useful  purposes  such  as  acceleration  or  deceleration  of 
flow,  prevention  of  separation,  and  the  like. 

Since  those  pioneering  years,  the  scientific/engineering  discipline  in  question  has 
been  going  through  a  continuous  process  of  maturation.  This  process,  however, 
has  been  continuously  and  systematically  marked  in  time  by  researchers’  com¬ 
plaints  about  the  unsatisfactory  state-of-the-art  of  the  theory.  Indeed,  notwithstand¬ 
ing  many  efforts,  and  the  voluminous  literature  generated  by  them,  to  confer  the 
discipline  the  status  of  being  firmly  established  on  physically  rigorous  and  consis¬ 
tent  foundations  freed  from  ad  hoc  assumptions,  progress  to  achieve  convergence 
to  that  goal  appears  today  not  completed  yet. 

The  study  presented  here  was  carried  out  in  the  context  of  a  research  activity  mo¬ 
tivated  by  renewed  interest  in  investigating  the  influence  that  electric  and/or  mag¬ 
netic  fields  can  exert  on  the  thermal  loads  imposed  on  a  body  invested  by  a  hyper¬ 
sonic  flow  [18,27-34].  In  this  regard,  spacecraft  thermal  protection  during  planetary 
(re)entry  represents  the  driving  engineering  application.  The  contents  of  the  study 
should  be  considered,  to  a  certain  extent,  a  systematic  reexamination  of  past  work 
complemented  with  somewhat  innovative  ideas.  The  aim  concentrates  on  the  for¬ 
mulation  of  a  consistent  set  of  governing  equations  in  open  form  apt  to  describe  the 
physical  phenomenology  comprising  the  hypersonic  flow  field  of  an  ionized  gas 
mixture  and  the  presence  of  the  electromagnetic  field.  The  discourse  opens  with 
stoichiometric  considerations  that  are  important  to  comprehend  how  specific  para¬ 
meters  of  electromagnetic  nature,  namely  electric-charge  density  and  conduction- 
current  density,  can  be  expressed  in  terms  of  variables  of  fluid-dynamics  nature. 
Subsequently,  the  governing  equations  of  the  flow  field  and  those  of  the  elec¬ 
tromagnetic  field  are  revisited  in  sequence;  differences  or  similarities  with  past 
treatments  are  pointed  out  and  discussed.  The  equations  governing  the  flow  field 
hinge  on  the  customary  balance  of  masses,  momenta  and  energies.  The  equations 
governing  the  electromagnetic  field  are  introduced  both  directly  in  terms  of  the 
Maxwell  equations  and  by  recourse  to  the  scalar  and  vector  potentials.  In  the  latter 
case,  the  convenience  of  adopting  the  Lorentz  gauge,  rather  than  the  magnetosta- 
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tic  gauge,  in  order  to  obtain  field  equations  with  favorable  mathematical  symmetry 
is  adequately  pointed  out.  Features,  limitations  and  approximations  implied  in  the 
open-form  governing  equations  are  explicitly  addressed.  Thermodynamics  aspects 
associated  with  the  necessity  to  assign  the  thermodynamic  model  of  the  gas  mix¬ 
ture  are  described  and  discussed.  The  theory  of  linear  irreversible  thermodynamics 
[20,21,23,35-37]  based  on  the  entropy-balance  equation  is  examined  for  the  pur¬ 
pose  of  obtaining,  consistently  with  the  presence  of  the  electromagnetic  field,  the 
phenomenological  relations  required  to  bring  into  a  mathematically  closed  form  the 
governing  equations.  Old  problems,  such  as  the  influence  of  medium  compressibil¬ 
ity  on  chemical-reaction  rates  or  the  importance  of  cross  effects  among  generalized 
fluxes  and  forces,  are  re-discussed;  additional  problems,  such  as  the  necessity  to 
take  into  account  the  tensorial  nature  of  the  transport  properties  because  of  the 
anisotropy  introduced  by  the  magnetic  field,  are  pointed  out.  A  non-conventional 
choice  of  first-tensorial-order  generalized  forces  and  corresponding  fluxes  is  pro¬ 
posed  which  appears  to  offer  more  simplicity  and  better  convenience  from  a  con¬ 
ceptual  point  of  view  when  compared  to  alternative  definitions  customarily  used  in 
the  literature. 

Polarization  and  magnetization  have  not  been  considered  in  this  study.  Setting  aside 
their  expected  negligibility  in  hypersonic  flows,  there  is  an  important  reason  behind 
that  choice.  The  inclusion  of  polarization  and  magnetization  effects  in  the  Maxwell 
equations  is  conceptually  (almost)  straightforward.  That,  however,  would  constitute 
only  a  unilateral  approach  to  the  physical  phenomenology.  Indeed,  the  important 
fact  should  not  be  overlooked  on  the  fluid-dynamics  side  that  not  only  body-force 
distributions  but  also  torque  distributions  exist  [21,38]  within  a  polarized  and  mag¬ 
netized  medium  subjected  to  the  action  of  the  electromagnetic  field.  Under  these 
circumstances,  the  velocity  vector  is  not  the  sole  kinematic  unknown  that  charac¬ 
terizes  the  flow  field;  the  specific  angular  momentum  [21,23,38]  of  matter  may  not 
identically  vanish  throughout  the  flow  field,  as  it  usually  happens  in  the  absence  of 
polarization  and  magnetization,  and  must  necessarily  be  taken  into  account  as  an 
additional  kinematic  unknown.  The  appearance  of  the  corresponding  balance  equa¬ 
tion  in  the  set  of  the  governing  equations  is  inescapable.  A  non-vanishing  specific 
angular  momentum  in  matter  can  have  far  reaching  consequences.  For  example, 
the  stress  tensor  loses  its  symmetry;  its  antisymmetric  part,  conjointly  with  the  an¬ 
tisymmetric  part  of  the  velocity-vector  gradient,  contributes  to  produce  entropy  and 
the  familiar  Newton  law  does  not  suffice  anymore  to  characterize  the  tensional  state 
in  the  medium.  Implications  of  energetic  nature  should  also  be  expected  because 
there  is  energy  associated  with  specific  angular  momentum;  in  addition,  the  polar¬ 
ization  and  magnetization  vectors  belong  to  the  set  of  the  thermodynamic  indepen¬ 
dent  state  parameters  [19,20].  These  and  similar  aspects  cannot  be  ignored  at  the 
moment  of  constructing  a  physically  rigorous  theory,  even  if  the  mentioned  effects 
may  turn  out  to  be  negligible  under  specific  flow  circumstances.  The  complexity 
of  the  physical  phenomenology  in  the  presence  of  polarization  and  magnetization 
increases  enormously  and  its  study  presupposes  a  degree  of  difficulty  which  can  be 
adequately  tackled  only  after  that  acquisition  of  solid  understanding  of  the  coupling 
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between  fluid  dynamics  and  pure  electromagnetic  field  has  been  secured.  The  latter 
constitutes  the  main  target  of  the  present  study  and  the  motivation  to  postpone  to 
future  investigations  the  behaviour  of  polarized  and  magnetized  media. 


2  Stoichiometric  aspects 


The  chemical  constituents  that  compose  an  ionized  gas  mixture  can  be  subdivided 
in  neutral  components,  ionized  components  and  free  electrons.  The  knowledge  of 
how  many  and  which  components  intervene  within  a  given  flow  problem  relies 
on  experimental  evidence  complemented  with  the  judicious  choice  dictated  by  the 
researcher’s  expertise.  The  acquisition  of  such  knowledge  is  sometimes  straightfor¬ 
ward,  sometimes  rather  involved;  in  any  case,  it  constitutes  a  problem  of  its  own. 
When  the  n  components  are  identified  then  it  is  possible  to  recognize  the  s  (s  <  n) 
reference  elements  that  participate  in  their  formation.  There  is,  obviously,  a  cer¬ 
tain  arbitrariness  in  the  qualification  of  the  reference  elements.  For  example,  either 
the  molecule  N2  can  be  considered  formed  by  putting  together  two  N  atoms  or, 
vice  versa,  the  atom  N  can  be  considered  formed  by  breaking  the  N2  molecule;  the 
role  of  reference  element  is  played  by  N  or  N2  in  the  former  or  latter  case,  respec¬ 
tively.  For  reasons  of  convenience,  the  reference  elements  are  enumerated  in  such 
a  way  that  the  first  (s  —  1)  are  the  true  atoms/molecules  E1,  E2, ...,  f^-1  and  the 
last  one  Es  is  the  electron  ‘e’  responsible  for  building  the  electric  charge  carried 
by  the  ionized  components  (if  any).  The  formation  concept  is  formalized  in  the 
chemical  formula  E^E^  ■  ■  ■  Esc.s  of  the  generic  component.  The  coefficient  oi:i 
represents  the  number  of  EJ  atoms/molecules  required  to  form  the  ?’-th  chemical 
component;  if  ar]  =  0  then  the  j-th  element  does  not  intervene  in  the  formation  of 
the  i-th  component  and  the  corresponding  symbol  (E])  is  dropped  from  the  chemi¬ 
cal  formula.  The  non- vanishing  coefficients  (j  =  1 .....  s  —  1 )  are  necessarily 
positive,  even  integers  if  the  elements  are  monatomic.  For  a  neutral  component,  the 
coefficient  als  is  identically  zero.  For  an  ionized  component,  the  coefficient  cris  is 
positive  or  negative  for  exceeding  or  missing  electrons  and  its  opposite  gives  the 
electric  charge  carried  by  the  component  molecule  as  an  integer  multiple  of  the 
electronic  charge.  In  a  more  common  formalism,  E*  is  replaced  by  superscripting 
the  chemical  formula  with  ‘+’  signs,  if  als  <  0,  or  ’  signs,  if  als  >  0,  in  number 
equal  to  als  .  The  coefficients  cr?J  can  be  grouped  together  to  compose  the  (n  x  s) 
formation  matrix.  The  rightmost  column  (j  =  s)  of  the  matrix  is  the  electric-charge 
column  which  appears  and  acquires  significance  exclusively  when  ionized  compo¬ 
nents  are  present  in  the  mixture.  The  coefficients  a,:i  permit  to  express  the  molar 
masses  Mj  of  the  components  as  linear  combinations 

Mi  —  ^2  aijM*  (1) 

3= 1 
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of  the  molar  masses  M*  of  the  elements  (M*  =  5.48579911  •  10'  '  kg  is  the  mass 
of  one  mole  of  electrons)  intervening  in  their  formation.  The  molar  electric  charges 
Qi  require  only  the  coefficients  in  the  electric-charge  column 


Qi  a 


(2) 


in  combination  with  the  electronic  charge  e  =  1.602176462  •  10”19C  and  the  Avo- 
gadro  number  NA  =  6.02214199  •  1023.  An  explicit  example  of  formation  matrix 
relative  to  a  seven-component  high-temperature  air  mixture  reads 


N 

0 

e 

N 

1 

0 

0 

0 

0 

1 

0 

e_ 

0 

0 

1 

NO 

1 

1 

0 

n2 
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0 

0 

02 

0 

2 

0 

NO+ 

1 

1 

-1 

n2 

o2 

e 

N 

1/2 

0 

0 

0 

0 

1/2 

0 

e“ 

0 

0 

1 

NO 

1/2 

1/2 

0 

n2 

1 

0 

0 

02 

0 

1 

0 

NO+ 

1/2 

1/2 

-1 

depending  whether  atoms  or  molecules  are  chosen  as  reference  elements.  In  this 
case,  there  are  n  =  7  components  formed  by  s  =  3  elements. 

Another  important  construct  is  the  stoichiometric  matrix  connected  with  the  r  chem¬ 
ical  reactions 

n  n 

£ ^2  [CTli  -  [CFli  k  =  l,2,...,r  (3) 

i=l  i=l 


that  can  occur  in  the  gas  mixture;  vkl ,  are  the  stoichiometric  coefficients  of  reac¬ 
tants  and  products,  respectively.  The  processes  formalized  in  Eq.  (3)  are  subjected 
to  component-mass  conservation 

n  n 

E  =  Y.  Mt  (4) 

i= 1  z=l 


After  defining  the  global  stoichiometric  coefficients  vkl  =  —  z/jj’,  Eq.  (4)  can  be 

recast  into  the  form 

n 

E^M*  =  0  (5) 

i=  1 
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The  coefficients  are  integer  numbers  and  are  conveniently  assembled  in  a  (r  x 
ri)  stoichiometric  matrix.  The  combination  of  the  formation  concept  embodied  in 
Eq.  (1)  and  the  component-mass  conservation  enforced  by  Eq.  (5)  allows  to  ob¬ 
tain  important  conditions  to  which  formation  and  stoichiometric  matrices  are  sub¬ 
jected  and  that  express  physically  the  element-mass  conservation.  The  substitution 
of  Eq.  (1)  into  Eq.  (5)  and  the  permutation  of  the  sum  operators  yields 

s  n 

M*  "kiCTij  =  o  (6) 

j= 1  i= 1 

Given  the  (mathematical)  arbitrariness  of  the  molar  masses  M*,  the  solution 

n 

uki°ij  =  0  (7) 

2=1 

is  the  sole  possibility  left  to  have  Eq.  (6)  identically  satisfied.  It  is  interesting  to 
notice  that  when  j  =  s,  and  taking  in  account  Eq.  (2),  Eq.  (7)  yields  the  electric- 
charge  conservation 

n 

vkiQi  =  o  (8) 

2=1 

across  the  given  chemical  reaction.  The  conservation  of  the  electric  charge  is,  there¬ 
fore,  not  an  independent  statement  but  follows  from  the  mass  conservation  relative 
to  the  electron  as  reference  element. 

The  formation  matrix  permits  to  express  composition  parameters,  and  their  proper¬ 
ties,  related  to  the  elements  in  terms  of  those  related  to  the  components.  The  basic 
relation,  in  this  regard,  is  the  one  that  connects  particle  number  of  the  elements  with 
particle  number  of  the  components 

n 

a;  (9) 

2=1 


From  Eq.  (9),  for  example,  one  obtains  similar  expressions  for  mass  fractions 


aj 


GO) 


and  partial  densities 


Pi 


(ID 


RTO-EN-AVT-162 


1  - 11 


Hypersonic-Flow  Governing  Equations  with  Electromagnetic  Fields 


ORGANIZATION 


The  electric  charge  per  unit  mass  also  follows  from  Eq.  (9)  after  setting  j  —  s  and 
multiplying  by  — e;  it  reads 


q  =  —  eNA-^j-  =  —eNA  Y" 

Ms  Mi 


(12) 


In  turn,  multiplication  of  Eq.  (12)  by  the  total-mass  density  p  provides  the  expres¬ 
sion  for  the  electric-charge  density 


Pc  =  pq 


-eNA 


P*s 

M* 


-e«AE 

i= 1 


pi_ 

Mi 


<7i. 


(13) 


Equation  (13)  is  an  important  relation.  It  establishes  a  first  necessary  link  between 
the  electromagnetic  side  (pc)  of  the  physical  phenomenology  and  its  fluid-dynamics 
counterpart  (p*  or  all  pt).  It  also  endorses  the  idea  that  the  electric-charge  density 
should  not  be  looked  at  as  a  basic  field  unknown  because  it  can  be  straightforwardly 
calculated  when  the  gas  mixture  composition  has  been  determined. 


3  Physical  significance  of  the  balance  equations 


In  view  of  the  analysis  in  the  following  sections,  it  appears  appropriate  to  dwell 
preliminarily  upon  an  important  aspect  related  to  the  physical  significance  of  the 
balance  equations  which  becomes  manifest  when  the  presence  of  the  electromag¬ 
netic  field  has  to  be  considered. 

It  is  a  recurrent  occurrence  in  the  mechanics  of  continuous  media  that  important 
equations  governing  the  dynamic  evolution  of  a  system,  namely  the  portion  of  the 
medium  contained  in  a  specified  control  volume,  are  developed  from  the  idea  of 
balancing  the  variations  of  the  extensive  properties  (mass,  momentum,  energy,  etc) 
that  characterize  the  macroscopic  state  of  the  system.  If  G  is  any  generic  extensive 
variable  owned  by  the  system  and  gv,  <bG,  gv  are  respectively  its  density,  flux  and 
production,  then  the  typical  balance  equation 

^  =  _V-S0  +  <7„  (14) 

is  the  translation  in  mathematical  language  of  the  basic  principle  [35,36,39]  affirm¬ 
ing  that  the  variable  G  can  vary  in  time  t  only  for  two  specific  reasons:  a)  an  ex¬ 
change  with  the  external  environment  and  b)  an  internal  production.  Equation  (14) 
constitutes  the  local  formulation  of  such  a  principle  and  establishes  a  formal  link 
between  time  variation  and  reasons  of  change.  Density  and  production  in  Eq.  (14) 
carry  the  subscript  v  to  emphasize  that  they  are  referred  to  unit  volume.  Feynman 
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provided  a  very  interesting  disquisition  in  his  famous  lectures  [40]  concerning  the 
physics  behind  Eq.  (14);  his  reasoning,  although  expounded  in  didactic  style,  is 
certainly  one  of  the  most  lucid  accounts  the  present  author  ever  had  the  opportunity 
to  read. 

The  aspect  related  to  Eq.  (14)  that  is  meant  to  be  pointed  out  here  regards  the  possi¬ 
bility  for  the  quantities  gv,  fly,  and  gv  of  being  attached  directly  to  space  instead  of 
being  associated  with  the  matter  occupying  that  same  space.  Such  an  occurrence  is 
somewhat  forgotten  in  traditional  fluid  dynamics  because  the  physical  variables  are 
all  associated  with  matter  in  that  case.  Indeed,  customary  practice  proceeds  one  step 
further  from  Eq.  (14)  by  introducing  density  g  =  gvj p  and  production  g  =  gv/p 
referred  to  unit  mass  and  by  separating  the  flux 


$g  =  pyg  +  Jg 


(15) 


in  a  convective  part,  associated  with  the  flow  velocity  v,  that  takes  care  of  the  trans¬ 
port  associated  with  the  macroscopic  motion  of  matter  and  a  diffusive  part  JG  which 
takes  care  of  everything  else.  Accordingly,  Eq.  (14)  becomes 


(16) 


Equation  (16)  is  the  stencil  that  embeds  all  the  governing  equations  belonging  to 
traditional  fluid  dynamics.  It  comes  to  no  surprise,  therefore,  that  the  mathematical 
structure  of  Eq.  (16)  has  stood  as  the  starting  point  in  computational  fluid  dynam¬ 
ics  (CFD)  from  which  all  efforts  towards  the  development  of  numerical  algorithms 
have  originated.  In  this  sense,  Eq.  (16)  has  undoubtedly  contributed  to  forging  the 
way  of  thinking  in  the  CFD  community.  Yet,  things  may  be  looked  at  from  a  dif¬ 
ferent  perspective  in  the  presence  of  the  electromagnetic  field.  Obviously,  the  ap¬ 
plicability  of  the  convection-diffusion  separation  [Eq.  (15)],  pertaining  to  the  flux 
fly  ,,  and  of  Eq.  (16)  still  survives  when  the  fluid-dynamics  field  and  the  electromag¬ 
netic  field  have  to  coexist.  However,  Eq.  (14)  can  also  play  a  role  if  the  definition 
of  global  variables,  namely  momentum  and  energy,  are  adequately  generalized  in 
a  manner  that  relaxes  the  unnecessary  conceptual  habit  of  matter  association;  then 
the  ensuing  equations  become  statements  of  conservation  (gv  =  0)  and,  in  so  doing, 
they  assume  a  mathematical  structure  that,  in  principle,  may  favorably  lend  itself 
to  a  more  simplified  numerical  analysis. 

In  the  following  sections,  the  formal  balance-equation  concept  will  be  explicited 
in  relation  to  the  fundamental  physical  quantities  mass,  momentum  and  energy 
in  order  to  formulate  a  consistent  set  of  governing  equations.  Concerning  the  lat¬ 
ter  two  quantities,  the  programme  will  be  carried  out  in  a  comparative  fashion  by 
confronting  the  fluid-dynamics  habitual  approach  with  the  novel  perspective  just 
discussed  and  brought  to  surface  by  the  presence  of  the  electromagnetic  field. 
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4  Mass-balance  equations 


The  standard  equations  balancing  the  component  masses 


dpi 

dt 


+  V  •  0;V) 


—  V  '  Jm,  +  Z 

k= 1 


i  —  1,  •  •  • ,  n 


(17) 


are  available  for  the  determination  of  the  gas-mixture  composition.  The  component- 
mass  diffusive  fluxes  Jm;  and  the  chemical-reaction  rates  require  the  assignment 
of  phenomenological  relations  (Secs.  8  and  9).  Other  variables  with  same  require¬ 
ment  will  be  encountered  in  the  sequel;  they  should  be  viewed  as  windows  through 
which  models,  describing  the  physical  behaviour  of  the  medium,  manifest  their 
influence  on  the  open-form  governing  equations.  The  component-mass  diffusive 
fluxes  are  linked  to  the  corresponding  diffusion  velocities 

Jm,  =  /Wi  (18) 


and  are  subjected  to  the  condition 

n  n 

ZJ  mi  =  ^2p%yfi  =  0  (19) 

2=1  2=1 

Equation  (19)  enforces  the  physical  fact  that  total  mass  cannot  diffuse.  In  other 
words,  there  are  only  n  —  1  independent  diffusive  fluxes  or  diffusion  velocities. 
Taking  into  account  Eq.  (5),  Eq.  (19)  and  mass  additivity 

n 

P  =  Epi  (20) 

2=1 

the  continuity  equation 

dp 

ft+V-  (pv)  =  0  (21) 

follows  from  the  summation  of  Eq.  (17)  on  the  subscript  i.  There  are  two  options 
for  the  determination  of  the  n  +  1  unknowns  p, ,  p.  The  most  straightforward  way 
would  seem  to  consist  in  the  selection  of  Eqs.  (17)  and  (20)  because  the  use  a 
very  simple  algebraic  equation,  as  Eq.  (20)  is,  is  appealing,  of  course.  There  is, 
however,  a  risk  in  doing  so  because  potential  inconsistencies  carried  into  Eq.  (17) 
by  phenomenological  relations  for  component-mass  diffusive  fluxes  and  chemical- 
reaction  rates  would  produce  inaccurate  partial  densities  which,  in  turn,  would  pass 
on  their  inaccuracy  to  the  total-mass  density  via  Eq.  (20).  The  alternative  way  to 
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proceed  could  be  to  replace  Eq.  (20)  with  Eq.  (21).  In  this  manner,  the  effect  of 
the  previously  mentioned  inconsistencies  is  somewhat  contained  because  neither 
partial  densities  nor  phenomenological  relations  are  explicitly  required  in  Eq.  (21). 
As  a  matter  of  fact,  Eq.  (20)  could  be  used  aside,  once  the  unknowns  have 
been  obtained,  as  a  sort  of  error  verifier.  The  drawback  of  this  approach  consists  in 
the  necessity  to  solve  an  additional  differential  equation  [Eq.  (21)]. 

The  chemical-reaction  rates  are  known  to  be  numerically  stiff  properties  to  deal 
with.  It  is,  therefore,  desirable  to  make  them  appear  as  sparingly  as  possible  in  the 
governing  equations.  To  this  aim,  simplification  can  be  achieved  to  some  extent  if 
the  element-composition  parameters  are  brought  into  the  picture.  Taking  into  ac¬ 
count  the  definition  of  element  partial  densities  [Eq.  (11)],  the  balance  equations  of 
the  element  masses  are  obtained  by  multiplying  Eq.  (17)  by  oy/M* /Mi  and  sum¬ 
ming  on  the  subscript  i\  they  read 


(22) 


The  element-mass  diffusive  flux  on  the  right-hand  side  of  Eq.  (22)  turns  out  to  be 
expressed  in  terms  of  the  component-mass  diffusive  fluxes  as 


(23) 


The  production  term  is  absent  in  Eq.  (22)  because  Eq.  (7)  makes  it  vanish  iden¬ 
tically.  Thus,  the  element  masses  are  conservative:  they  cannot  be  either  created 
or  destroyed,  regardless  of  the  reactive  mechanisms  at  work  in  the  gas  mixture. 
This  occurrence  suggests  an  advantageous  manoeuvre  to  limit  the  appearance  of 
the  chemical-reaction  rates.  The  idea  is  to  relinquish  as  unknowns  the  last  s  partial 
densities  pi  and  replace  them  with  the  partial  densities  p*;  at  the  same  time,  the  last 
s  equations  of  the  set  (17)  are  replaced  with  the  set  (22).  In  this  way,  the  number  of 
differential  equations  is  unchanged  but  the  chemical-reaction  rates  appear  only  in 
n  —  s  equations.  The  s  relinquished  partial  densities  p*  can  be  expressed  in  terms 
of  the  first  n  —  s  partial  densities  pt  and  of  the  s  partial  densities  p*  from  Eq.  (11) 
after  expanding 


3  =  ■  ■  ■ ,  « 


(24) 


and  re-arranging 


n 


E 


3  = 


(25) 


i=n—s+ 1 
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Equation  (25)  represents  an  algebraic  system  of  s  equations  for  the  s  relinquished 
partial  densities  p*.  It  can  be  formally  solved  in  the  following  manner.  First,  intro¬ 
duce  for  brevity  the  (s  x  s)  coefficient  matrix 

M* 

aij  =  aijl\t  i  =  n-s+l,---,n;j  =  l,---,s  (26) 


and  the  (1  x  s)  known-term  array 

n~s  o 

p]  = $  -  E  ,  j  =  i> •  •  • » s 


2=1 


so  that  Eq.  (25)  can  be  recast  in  the  standard  form 

n 

El  Pi&ij  =  Pj  j  =  1)  ■  ■  ■  J  s 

i=n-s-\-l 


(27) 


(28) 


The  matrix  atJ  can  be  inverted  once  and  forever  when  the  formation  matrix  <Jij  and 
the  element  molar  masses  are  known.  Then,  the  formal  solution  of  Eq.  (28)  is 

S 

Pi  =  J2  P*j  <*ji  i  =  n-s  + 1,  — ,  n  (29) 

j= i 


Another  important  aspect  to  look  at  in  connection  with  the  elements  is  the  balancing 
of  the  electric  charge.  The  electric -charge  balance  equation  is  not  an  independent 
statement  but  is  embedded  in  Eq.  (22)  when  particularized  to  the  case  of  the  electron 
element.  Indeed,  setting  j  =  s  in  Eq.  (22)  and  multiplying  it  by  —eNA/M*  [see 
Eq.  (13)]  yields  the  fluid-dynamics  styled  equation 

do 

^  +  V-(Pcv)  =  -V-  JQ  (30) 


in  which,  taking  into  account  Eq.  (23)  with  j  =  s,  the  electric -charge  diffusive  flux 
turns  out  to  be  expressed  by  the  following  linear  combination 


Jq  = 


-eNJ 


2=1 


(31) 


of  the  component-mass  diffusive  fluxes  of  the  electrically  charged  components 
{a iS  7^  0).  It  is  important  to  notice  that  Eqs.  (30)  and  (31)  warn  against  any  pre¬ 
sumptive  imposition  of  charge  neutrality  (pc  =  0)  throughout  the  flow  field;  even 
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if  there  are  zones  in  which  the  gas  mixture  is  electrically  neutral,  the  mass  diffu¬ 
sion  of  ionized  components  and  free  electrons  works  towards  the  removal  of  such 
a  condition.  Rephrasing  Eq.  (30)  in  electromagnetic-theory  style 

dp 

-^  +  V-(pcv  +  JQ)  =  0  (32) 

leads  to  the  identification  of,  in  the  corresponding  parlance,  the  electric-current 
density 

j  =  Pc-V  +  JQ  (33) 

and  to  the  recognition  of  its  separability  in  convection-current  (pcv)  and  conduction- 
current  (JQ)  densities.  These  are  the  sole  contributions  that  need  to  be  accounted 
for  in  the  absence  of  polarization  and  magnetization.  Equation  (31)  represents  the 
other  important  relation  that  establishes  a  second,  and  final,  necessary  link  between 
electromagnetism  (JQ)  and  fluid  dynamics  (all  Jmi).  Both  Eq.  (13)  and  Eq.  (31) 
converge  into  the  definition  provided  by  Eq.  (33)  and,  in  so  doing,  enforce  the 
unambiguous  assertion  that  the  electric-current  density  is  specified  entirely  in  terms 
of  variables  of  fluid-dynamics  nature.  Also,  the  dependence  expressed  in  Eq.  (31) 
clearly  shows  that  it  is  not  necessary  to  pursue  an  independent  phenomenological 
relation  for  the  conduction-current  density  because  the  latter  descends  naturally 
from  the  knowledge  of  the  phenomenological  relations  for  the  component-mass 
diffusive  fluxes.  It  will  be  seen  in  Sec.  9.4  how  the  famous  Ohm  law  and  additional 
effects  of  thermodynamic  origin  arise  naturally  in  this  way. 


5  Electromagnetic-field  equations 


The  essence  of  electromagnetism  finds  its  deepest  representation  in  the  differential 
equations  that  govern  the  electromagnetic  field,  namely  the  well  known  Maxwell 
equations.  The  body  of  didactic  literature  on  this  subject  matter  is  enormous  and  the 
theory  can  very  well  be  considered  consolidated  on  solid  foundations.  The  contents 
of  this  section  take  advantage  mainly  from  Feynman’s  lectures  [40];  Maxwell’s  fun¬ 
damental  treatise  [41,42]  together  with  the  textbooks  written  by  Mpller  [43],  Lor- 
rain  and  Colson  [44],  Persico  [45],  Tolman  [46],  and  Pauli  [47]  were  also  helpful. 
Notwithstanding  the  satisfactory  state-of-the-art  of  the  theory,  there  is  one  peculiar 
aspect  of  electromagnetism  that  always  deserves  extreme  care  and  attention:  the 
choice  of  the  physical  units.  In  SI  units,  the  Maxwell  equations  read 


V  -E=- 

(34) 

£o 

V  B  =  0 

(35) 
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„  OB 

VXE=-^ 

(36) 

.  0E 

e0c  VxB=j  +£o~fa 

(37) 

The  electric-charge  and  electric-current  densities  represent  the  channels  through 
which  the  coupling  between  fluid-dynamics  field  and  electromagnetic  field  be¬ 
comes  manifest  [recall  Eq.  (13),  Eq.  (31)  and  Eq.  (33)].  The  constants  c  and  e0 
are  respectively  the  velocity  of  light  (299792458  m-s-1)  in  and  the  dielectric  con¬ 
stant  (8.854187817 •  10“12  F-m”1)  of  vacuum.  In  principle,  Eqs.  (36)  and  (37)  are  all 
that  is  required  to  associate  with  the  fluid-dynamics  equations  in  order  to  determine 
simultaneously  electric -field  intensity  E  and  magnetic  induction  B.  However,  their 
mathematical  structure  is  substantially  distinct  from  the  habitual  fluid-dynamics 
stencil  [Eq.  (16)].  A  widespread  practice  [34,48-55]  that  aims  to  derive  and  use  an 
equation  with  more  CFD-suitable  form  is  based  on  the  adaptation  of  Eq.  (36)  fol¬ 
lowing  the  neglect  of  the  displacement-current  density  (e0dE/dt)  in  Eq.  (37)  and 
the  assumed  validity  of  the  generalized  Ohm  law 

JQ  =  X,  (E  +  v  x  B)  (38) 


although  with  a  scalar  electrical  conductivity  Xc.  The  method  leads  to  an  algebraic 
relation  for  the  electric  field 

2 

E=— VxB--v-vxB  (39) 

Xc  Xc 

and  to  the  so-called  magnetic-induction  equation 

f  c^1  f  C‘ ^ 

—  +  V  ■  (vB)  =  V  •  (Bv)  +  — V2B  +  ^VXe  x  (V  x  B) 

Ot  Xc  Xc 

+  ^Vxv-vxV  (40) 

The  details  of  the  derivation  of  Eq.  (40)  are  given  in  appendix  A.  Further  simpli¬ 
fied  forms  in  the  event  of  electric-charge  neutrality  (pc  ~  0)  or  uniform  electrical 
conductivity  (Xe  ~const)  are  easily  deduced.  Equation  (40)  looks  certainly  attrac¬ 
tive  from  a  numerical  point  of  view  because  its  structure  reflects  perfectly  that  of 
Eq.  (16).  In  this  way,  the  solution  of  the  electromagnetic  field  is  brought  within  the 
reach  of  familiar  algorithms  in  CFD.  At  the  same  time,  the  idea  of  magnetic-field 
convection  is  favoured  to  find  its  way  into  the  picture  of  the  physical  phenomenol¬ 
ogy.  The  computational  fluid  dynamicist  is  most  likely  satisfied  with  this  situation 
because  he  is  provided  with  an  additional  instrument  [Eq.  (40)]  which  features  the 
same  mathematical  characteristics  of  familiar  tools,  namely  the  flow-field  equa¬ 
tions  without  the  electromagnetic  field.  He  can,  then,  proceed  to  calculate.  Seen 
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from  the  perspective  of  the  hard  efforts  and  time  invested  in  the  development  of 
numerical  schemes,  this  attitude  is  comprehensible.  Yet,  the  theoretical  fluid  dy- 
namicist  would  feel  concerned  about  the  same  situation  because  Eq.  (40)  and  its 
associated  interpretation  are  very  much  in  contrast  with  the  physical  fact  that  the 
electromagnetic  field  is  attached  to  space  regardless  of  the  matter  flowing  through 
that  same  space.  In  this  regard,  he  may  ponder  about  the  physical  significance  of 
Eq.  (40),  asking  important  questions  such  as:  can  the  vector  B  be  interpreted  as  the 
volume  density  of  some  extensive  property  of  the  matter  moving  in  the  space  oc¬ 
cupied  by  the  electromagnetic  field?  If  there  is  such  a  property  then  can  the  diadic 
tensor  Bv  and  the  remaining  terms  on  the  right-hand  side  of  Eq.  (40)  be  interpreted 
as,  respectively,  its  diffusive  flux  and  production?  The  difficulty  in  finding  con¬ 
vincing  answers  suggests  a  critical  scrutiny  of  the  assumptions  on  which  Eq.  (40) 
is  built.  The  neglect  of  the  displacement-current  density  is  justifiable  in  circum¬ 
stances  of  not  rapidly  varying  electric  field  but  it  is  still  an  undesirable  limitation  at 
the  moment  of  constructing  a  general  theory.  The  generalized  Ohm  law  [Eq.  (38)] 
is  more  prone  to  criticism.  Concern  about  its  applicability  is  not  a  novelty  and  was 
explicitly  raised  long  time  ago  by  Maxwell  [41,42]  and  emphasized  in  more  recent 
times  by  Napolitano  [11,16],  Pai  [18]  and  Sedov  [21].  The  major  hurdle  to  accept 
is  the  fact  that  the  applicability  of  Eq.  (40),  which  is  a  governing  equation,  is  subju¬ 
gated  to  the  validity  of  Eq.  (38),  which  is  a  phenomenological  relation.  This  levies 
a  serious  toll  on  the  generality  of  the  ensuing  theory  because  the  latter  becomes 
medium-dependent.  Equation  (17),  for  example,  is  medium-independent  because  it 
remains  applicable  regardless  of  the  phenomenological  relations  assumed  for  the 
component-mass  diffusive  fluxes  and  chemical-reaction  rates.  This  is  not  the  case 
for  Eq.  (40).  What  happens  if  the  tensorial  nature  of  the  electrical  conductivity,  a 
feature  already  discussed  by  Maxwell  in  1873,  cannot  be  neglected  or,  worst,  if  the 
medium  does  not  comply  with  Eq.  (38)?  Indeed,  and  just  to  mention  an  example, 
Ohm  law  [Eq.  (38)]  becomes  meaningless  for  a  polarizable  and  magnetizable  neu¬ 
tral  gas.  There  can  be  no  electrical  conduction  (J0  =  0)  in  such  a  gas  because  free 
electric  charges  are  absent;  yet  there  is  an  electrical-current  density 

<9P 

3  =  ft  +  V  x  M  (41) 

produced  by  the  polarization  P  and  magnetization  M  of  the  gas.  In  this  case,  the 
whole  edifice  built  on  the  magnetic-induction  equation  [Eq.  (40)]  must  be  thrown 
away  because  absolutely  inapplicable  and  a  new  theory  must  be  constructed  afresh. 
These  arguments  may  appear  irrelevant  to  the  computational  fluid  dynamicist  who 
is  interested  mainly  in  numerical  algorithms  but  for  the  theoretical  fluid  dynamicist 
they  are  strong  reasons  of  concern  that  originate  from  the  awareness  of  operating 
on  the  basis  of  a  theory  whose  solidity  may  be  compromised  in  unforeseeable  and 
uncontrollable  particular  situations. 

Another  exploitable  method  takes  advantage  of  the  scalar  and  vector  potentials 
0,  A  often  used  in  electromagnetism.  The  knowledge  of  the  potentials  implies  that 
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of  electric-field  intensity  and  magnetic  induction  because  the  latter  vectors  follow 
from  the  definitions 


1 

<1 

-e- 

1 

<-+ 1  >- 

(42) 

V  x  A 

(43) 

The  substitution  of  Eqs.  (42)  and  (43)  into  Eqs.  (36)  and  (37)  leads  to  the  following 
nicely  symmetrical  field  equations 


1  <920 

c2  dt 2 

1  d2A 
c2  dt2 


V20+^ 

£0 


V2A  + 


(44) 

(45) 


The  details  of  the  derivation  are  provided  in  appendix  B.  The  mathematical  sym¬ 
metry  of  Eqs.  (44)  and  (45)  is  strongly  dependent  on  the  adoption  of  the  condition 


1  <90 

c2  dt 


+  V  •  A  =  0 


(46) 


known  as  Lorentz  gauge  [40,47].  A  similar  approach  was  already  pursued  by  Burg¬ 
ers  [14]  and  Pai  [18]  who,  however,  opted  for  the  typical  magnetostatic  gauge 

V  ■  A  =  0  (47) 


Instead  of  Eqs.  (44)  and  (45),  they  obtained  two  much  more  complicated  highly 
cross-coupled  field  equations  in  which  terms  involving  0  and  A  appear  simultane¬ 
ously  in  both  equations. 

Equations  (44)  and  (45)  indicate  explicitly  the  wave-like  evolution  taking  place  in 
the  electromagnetic  field  and  how  that  is  influenced  by  the  presence  of  matter 
through  the  electric-charge  and  electric-current  densities.  The  equations  reduce  to 
the  Poisson  equation  in  steady- state  circumstances.  It  may  be  asked  what  is  the 
gain  of  using  Eqs.  (44)  and  (45)  rather  than  Eqs.  (36)  and  (37)  or  Eq.  (40).  First 
of  all,  Eqs.  (44)  and  (45)  are  four  scalar  differential  equations  instead  of  the  six 
represented  by  Eqs.  (36)  and  (37).  Moreover,  they  are  general  and  independent  of 
the  medium  in  so  far  as  they  are  unaffected  by  arguments  related  to  importance  or 
disregard  of  the  displacement-current  density  and  as  they  need  no  appeal  to  any 
phenomenological  relation  to  provide  reason  for  their  existence.  It  is  true  that,  once 
again,  Eqs.  (44)  and  (45)  do  not  reflect  the  structure  of  Eq.  (16)  and,  therefore, 
they  presuppose  the  necessity  to  develop  new  numerical  algorithms  for  their  simul¬ 
taneous  solution  with  the  fluid-dynamics  equations.  On  the  other  hand,  they  are 
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equations  of  the  mathematical  physics  which  have  been  studied  numerically  since 
long  time  and  for  whose  properties  a  huge  body  of  knowledge  and  understanding 
has  been  accumulated. 


6  Momentum-balance  equations 


The  determination  of  the  velocity-vector  field  occurs  via  the  equation  balancing  the 
momentum  associated  with  the  matter  flowing  in  the  control  volume.  For  reasons 
that  will  appear  evident  soon,  it  is  appropriate  to  emphasize  the  association  to  mat¬ 
ter  by  systematically  referring  to  this  quantity  with  the  term  matter  momentum.  In 
the  presence  of  the  electromagnetic  field,  its  balance  equation  assumes  the  form 

^  +  V  •  (pv  v)  =  V  ■  T  +  pcE  +  j  x  B  (48) 

at 

As  in  traditional  fluid  dynamics,  the  matter-momentum  diffusion  is  characterized 
by  the  stress  tensor  T  which  requires  the  assignment  of  a  phenomenological  re¬ 
lation  and,  as  anticipated  in  Sec.  1,  preserves  its  feature  of  being  a  symmetrical 
tensor  in  the  absence  of  polarization  and  magnetization.  The  gravitational  contri¬ 
bution  to  the  body  force  on  the  right-hand  side  of  Eq.  (48)  has  been  omitted  for 
compatibility  with  the  typical  circumstances  settling  in  in  hypersonic  regime  that 
presuppose  the  negligibility  of  gravitational  effects  with  respect  to  those  due  to  the 
tensional  state  of  the  medium.  As  a  matter  of  fact,  the  inclusion  in  the  discourse  of 
a  (Newtonian)  gravitational  field  is  conceptually  straightforward  because  the  grav¬ 
itational  body-force  term  can  be  treated  similarly  to  the  electric  counterpart  (pcE) 
and  made  fit  smoothly  in  the  equation  framework  described  in  the  sequel.  On  the 
other  hand,  the  emphasis  of  the  present  context  addresses  the  importance  of  the 
electromagnetic  field;  the  presence  of  a  gravitational  field  would  only  burden  the 
equations  with  unnecessary  additional  terms  whose  inclusion  would  not  change  at 
all  the  considerations  that  will  follow  and  the  conclusions  that  will  be  drawn. 

The  electromagnetic  field  produces  matter  momentum  through  the  body  force 

g„  =  P.E  +  j  x  B  (49) 


This  is  the  term  to  which  the  statement  of  Resler  and  Sears  [9]  quoted  in  the  in¬ 
troduction  refers  to  and  that  is  responsible  for  a  variety  of  new  effects  substantially 
inimaginable  in  traditional  fluid  dynamics.  Under  the  action  of  the  electromagnetic 
field,  for  example,  the  mass  diffusion  of  the  electrically  charged  components  takes 
an  active  role  in  affecting  the  motion  of  the  fluid  particles  because  it  enters  ex¬ 
plicitly  into  the  equation  of  their  motion  [Eq.  (48)]  via  the  electric-current  density 
[see  Eqs.  (3 1)  and  (33)] .  Without  the  electromagnetic  field,  mass  diffusion  produces 
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only  an  indirect  effect  on  the  dynamics  of  the  flow  field  through  a  thermodynamic 
pathway  that  involves  the  gas-mixture  composition  and,  subsequently,  the  pressure 
distribution;  the  latter,  in  turn,  represents  a  substantial  contribution  to  the  build-up 
of  the  stress  tensor. 

Equation  (48)  is  a  necessary  and  sufficient  equation  qualified  for  inclusion  in  the 
governing  set;  one  could  be  satisfied  with  its  availability.  Nevertheless,  there  are 
more  interesting  features  of  the  physical  phenomenology  that  await  to  be  unrav¬ 
elled.  Whether  it  may,  perhaps,  appear  a  somewhat  fortuitous  circumstance  or  it 
could  be  looked  at  as  the  manifestation  of  something  of  deep  physical  significance, 
it  is  certainly  interesting  that  the  Maxwell  equations  [Eqs.  (34)  and  (37)]  allow  a 
very  useful  transformation  [43,45-47]  of  the  electromagnetic  body  force  (49).  In  a 
single  stroke,  this  transformation  provides  evidence  of  the  existence  of  momentum 
associated  with  the  electromagnetic  field,  namely  the  electromagnetic  momentum, 
and  leads  to  the  formulation  of  its  balance  equation.  The  mathematical  details  are 
given  in  appendix  C.  The  final  outcome  from  the  mentioned  transformation  pro¬ 
vides  the  electromagnetic  body  force  in  the  form 

pcE  +  j  x  B  =  V  •  XM  —  ^(£oE  x  B)  (50) 


In  Eq.  (50),  the  tensor  TM  represents  the  following  combination 

Tm  =  £o(EE  -  l-E2U)  +£0C2(BB  -  X-B2 U)  (51) 

of  electric-field  intensity,  magnetic  induction  and  unit  tensor  U.  It  is,  therefore,  a 
symmetric  tensor.  Equation  (50)  can  be  simply  overturned  as 

^(£0E  x  B)  =  V  •  Tm  -  (pcE  +  j  x  B)  (52) 

to  match  exactly  the  structure  of  Eq.  (14).  Hence,  Eq.  (52)  is  a  balance  equation. 
It  shows  unequivocally  the  existence  of  electromagnetic  momentum  distributed  in 
space  with  density  e0E  x  B  and  transported  through  space  with  flux  — TM.  The  sym¬ 
metric  tensor  TM  plays  in  Eq.  (52)  the  same  role  fulfilled  by  the  stress  tensor  in 
Eq.  (48)  and,  for  this  reason,  it  is  suggestively  named  as  Maxwell  stress  tensor. 
Equation  (52)  highlights  in  an  evident  manner  also  that  the  transport  of  electro¬ 
magnetic  momentum  takes  place  through  space  exclusively  in  consequence  of  the 
presence  of  the  electromagnetic  field  and  bears  no  relation  whatsoever  with  the 
matter  transported  through  that  same  space.  It  is  important  to  keep  in  mind  that, 
although  very  useful,  Eq.  (52)  is  not  a  new  independent  equation.  In  principle,  it 
can  replace  one  of  Eqs.  (36)  and  (37)  but  it  does  not  say  anything  more  that  is  not 
already  contained  in  the  Maxwell  equations.  The  nice  features  of  Eq.  (52)  consist 
in  its  balance-equation  structure  and  that  it  fulfills  the  task  of  permitting  a  deep 
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insight  towards  the  understanding  of  the  existence  of  important  properties  of  the 
electromagnetic  field  that  are  not  immediately  recognizable  from  the  equations  of 
electromagnetism  as  given  by  Maxwell. 

Another  aspect  worth  of  attention  is  that  all  the  electromagnetic  momentum  that 
disappears  locally  turns  out  to  reappear  as  matter  momentum  or  viceversa.  This  is 
the  obvious  conclusion  ensuing  from  the  appearance  of  the  electromagnetic  body 
force  both  in  Eq.  (48)  and,  with  changed  sign,  in  Eq.  (52).  In  other  words,  the  sum  of 
the  two  forms  of  momentum  cannot  be  produced,  either  created  or  destroyed.  Thus, 
the  global  momentum  pv+e0E  x  B  is  a  conservative  property  of  the  physical  system 
composed  by  the  conjoint  fluid-dynamics  and  electromagnetic  fields.  This  profound 
characteristic  of  the  physical  phenomenology  is  brought  to  surface  by  summing 
together  Eqs.  (48)  and  (52)  to  obtain  the  balance  equation  of  total  momentum 

-^-(pv  +  e0E  x  B)  =  —V  ■  (pvv  —  T  —  Tm)  (53) 

Equation  (53)  is  equivalent  to  Eq.  (48)  and  constitutes  a  valid  and,  perhaps,  more 
convenient  alternative  at  the  moment  of  performing  numerical  calculations  because 
it  is  not  burdened  by  the  presence  of  any  production  term. 


7  Energy-balance  equations 

7.1  Preliminary  considerations 


The  prerequisite  steps  in  the  formulation  of  balance  equations  related  to  the  concept 
of  energy,  in  its  entirety  and  in  its  variety  of  kinds,  are  the  identification  of  the  forms 
that  play  a  role  within  a  specific  physical  phenomenology  and  the  recognition  of 
the  sum  of  those  forms  as  the  total  energy.  In  turn,  the  subduing  of  the  latter’s 
production  ev  per  unit  volume  and  time  to  the  famous  principle  of  conservation 
(ev  =  0)  leads  to  the  deduction  and,  at  the  same  time,  the  physical  interpretation 
of  interesting  and  important  features  related  to  the  possible  mechanisms  of  energy 
conversion. 

In  the  absence  of  electromagnetic  fields,  the  typical  forms  of  energy  that  intervene 
in  hypersonic  regime  are  the  kinetic  energy  possessed  by  the  fluid  particles  as  a 
consequence  of  their  macroscopic  motion  and  the  internal  energies  distributed  over 
the  molecular  degrees  of  freedom  of  the  components.  Energy  ascribed  to  intermole  - 
cular  interactions  is  systematically  neglected.  Gravitational  energy  is  not  admitted 
in  the  picture  for  the  same  reason  of  negligibility  adduced  in  Sec.  6  to  justify  the 
omission  of  the  gravitational-field  contribution  to  the  body  force  on  the  right-hand 
side  of  Eq.  (48).  Under  the  assumed  circumstances,  the  sum  of  kinetic  energy  and 
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internal  energies  constitutes  the  total  energy  and,  as  such,  that  sum  acquires  the 
prerogative  of  being  conservative.  The  situation  changes  drastically  in  the  presence 
of  the  electromagnetic  field.  The  kinetic  energy  and  the  internal  energies  associated 
with  matter  are  still  part  of  the  scene  but  their  sum,  which  will  be  referred  to  as  mat¬ 
ter  energy  for  consistency  with  the  terminology  introduced  in  Sec.  6  when  dealing 
with  momentum,  does  not  exhaust  the  totality  of  forms.  The  electromagnetic  field 
possesses  energy  in  the  same  way  as  it  does  for  momentum.  It  turns  out,  therefore, 
that  there  is  a  further  kind  of  energy  to  account  for:  the  electromagnetic  energy.  It  is 
the  sum  of  matter  energy  and  electromagnetic  energy  to  provide  the  total  energy  in 
this  case  and  to  be  characterized  by  a  vanishing  production.  It  will  be  shown  in  the 
sequel  that,  once  again,  the  Maxwell  equations  [Eqs.  (36)  and  (37)]  and  their  ade¬ 
quately  manipulated  blend  with  the  balance  equations  of  kinetic  energy  and  internal 
energies  play  a  fundamental  role  in  the  achievement  of  the  outlined  understanding 
of  the  physical  situation. 


7.2  Kinetic  energy 


The  kinetic-energy  balance  equation 

+  V  '  G°yv)  =  V  •  (T  •  v)  -  T  :  Vv  +  pcv  •  E  -  JQ  •  v  x  B  (54) 

descends  straightforwardly  from  that  of  matter  momentum  simply  by  scalar-multi- 
plying  both  sides  of  Eq.  (48)  by  the  velocity  vector  and  by  rearranging  the  resulting 
right-hand  side  to  reflect  the  structure  of  Eq.  (16).  Inspection  of  Eq.  (54)  indicates 
at  once  kinetic-energy  diffusive  flux  and  production.  The  latter  comprises  the  ha¬ 
bitual  contribution  that  includes  the  combined  action  of  medium  deformation  and 
tensional  state,  and  a  contribution  originating  from  the  existence  of  the  electromag¬ 
netic  field.  With  regard  to  this  additional  contribution,  the  magnetic  part  contains 
only  the  conduction-current  density.  This  is  the  obvious  consequence  of  the  or¬ 
thogonality  [v  ■  (pcv  x  B  )  =  0]  between  the  velocity  vector  and  the  part  of  the  body 
force  in  Eq.  (49)  containing  the  convection-current  density  that  appears  explicitly 
after  expanding  the  electric-current  density  according  to  Eq.  (33).  Equation  (54)  is, 
clearly,  not  an  independent  equation;  it  merely  represents  the  projection  of  Eq.  (48) 
along  the  local  direction  of  the  instantaneous  streamlines  of  the  flow  field. 


7.3  Internal  energy 


From  a  thermodynamic  point  of  view,  the  ionized  gas  mixture  of  interest  in  the 
present  context  has  to  be  considered  as  a  composite  system  whose  subsystems,  rep¬ 
resented  by  the  molecular  degrees  of  freedom  possessed  by  the  components,  are  in 
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disequilibrium  with  respect  to  mass  exchanges  (chemical  reactions)  and  energy  ex¬ 
changes  (thermal  relaxations)  [39].  In  other  words,  the  internal  energies  distributed 
over  the  molecular  degrees  of  each  component  must  be  introduced  and  accounted 
for  separately.  In  general,  the  i-th  component  owns  independent  molecular  de¬ 
grees  of  freedom  (5  —  1,  •  •  • ,  ii)  and  the  4-th  degree  of  freedom  features  its  private 
internal  energy  Uis.  It  appears  worthwhile,  incidentally,  to  mention  that  the  prob¬ 
lem  of  the  explicit  separation  of  the  molecular  degrees  of  freedom  in  independent 
entities  at  the  level  of  the  internal  Schrodinger  equation  of  the  molecules  is  still  an 
open  issue  in  demand  of  satisfactory  resolution  and  is  systematically  glossed  over 
by  making  recourse  to  the  poor,  and  incorrect,  classical  separation  in  electronic, 
vibrational,  rotational  (and  etc)  molecular  degrees  of  freedom.  This  is  certainly 
a  gap  that  calls  for  enhancement  of  basic  understanding  achievable  only  through 
advanced  research.  The  pointed-out  limitation,  however,  does  not  prevent  the  de¬ 
velopment  of  a  formal  equation  framework.  The  independence  of  the  degrees  of 
freedom  implies  the  additivity  of  the  internal  energies 

n  li 

U  =  J2J2UiS  (55) 

i= 1  <5=1 


Equation  (55)  gives  the  internal  energy  of  the  gas  mixture  and  can  be  conveniently 
rephrased  in  terms  of  specific  quantities  as 

n  ii 

/W  =  EE  PiuiS  (56) 

i= 1  <5=1 


The  specificity  of  u^s  refers  to  the  unit  mass  of  the  /’-th  component  while  the  speci- 

n 

ficity  of  v  refers  to  the  unit  total  mass.  On  the  fluid-dynamics  side,  the  E  (,  specific 

i= 1 

internal  energies  are  unknowns  of  the  flow  field  and  their  determination  can  be 

n 

achieved  through  the  following  E  U  balance  equations 

Z=1 


dp^s 

dt 


+  V  •  {piUifiv) 


^  '  Jui5  T  Uv,iS 


(57) 


<5  =  1,  •  ■  •  i  =  1,  ■  ■  •  ,n 


The  diffusive  fluxes  and  productions  appearing  on  the  right-hand  side  of  Eq.  (57) 
require  the  assignment  of  phenomenological  relations.  With  regard  to  the  produc- 

n 

tions,  it  will  be  shown  in  Sec.  7.4  that  only  E  —  1  of  them  are  independent  in 

Z— 1 

consequence  of  the  principle  of  total-energy  conservation.  On  the  thermodynam¬ 
ics  side,  the  specific  internal  energies  Ui$  are  linked  to  the  Helmholtz  potentials 
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fis(Tis,  Vi)  that  describe  the  thermodynamic  behaviour  of  the  molecular  degrees  of 
freedom  via  their  dependence  on  the  corresponding  temperatures  and  on  the 
specific  volumes  of  the  components.  The  link  takes  the  form 


(58) 


5  =  1,  •  •  •  ,£i\i  =  1,  •  •  •  ,n 


n 


Each  of  the  ^  ll  thermodynamic  relations  (58)  provides  the  functional  dependence 


2=1 


to  obtain  the  temperature  Tt$  for  prescribed  specific  internal  energy  Ui$  and  specific 
volume  (vi  =  1  j  pi)  of  the  i-th  component.  The  set  of  the  functions  ./'„j  (T>rs ,  v,  )  char¬ 
acterizes  the  global  thermodynamic  model  [39]  of  the  gas  mixture.  Their  explicit 
determination  presupposes  the  knowledge  of  appropriate  partition  functions  [56,57] 
whose  construction,  in  turn,  belongs  to  the  domain  of  statistical  thermodynamics 
(Refs.  [56-59]  and  references  therein).  It  ought  to  be  remarked  that  the  described 
equation  scheme  is  founded  on  the  assumption  that  the  population  distributions 
over  the  quantum-energy  states  associated  with  the  molecular  degrees  of  freedom 
can  be  represented  in  analytical  form,  the  Boltzmann  distribution  being  a  particular 
case.  This  assumption  is  critical  for  the  effectiveness  of  Eq.  (57)  and  the  validity 
of  Eq.  (58).  Circumstances  cannot  be  excluded  in  which  this  assumption  becomes 
untenable.  In  that  case,  a  deeper  characterization  of  the  thermal  relaxations,  with 
repercussions  on  the  chemical  kinetics  of  the  gas  mixture,  becomes  necessary  be¬ 
cause  the  quantum-state  populations  are  themselves  unknowns  subjected  to  balance 
equations  that  deal  with  state-to-state  exchanges  of  energy  and  mass.  A  substantial 
body  of  works  (Refs.  [60-64]  and  references  therein)  addressing  the  state-to-state 
phenomenology  has  been  growing  recently  but  the  methods  elaborated  so  far  are 
not  yet  completely  free  from  difficulties  of  conceptual  and  computational  nature. 
However,  these  difficulties  notwithstanding,  experimental  and  computational  evi¬ 
dence  (Refs.  [65,66]  and  references  therein)  of  the  existence  of  non- analytical  dis¬ 
tributions  of  the  quantum-state  populations  points  towards  the  conclusion  that  the 
avenue  of  state-to-state  thermal  kinetics  certainly  deserves  to  be  explored  with  vig¬ 
orous  effort  for  reasons  of  both  scientific  and  engineering  interest.  This  topic  will 
not  be  elaborated  further  here  because  it  is  beyond  the  scope  of  the  present  context. 
Interested  readers  are  referred  to  the  cited  literature. 

Taking  into  account  the  additivity  [Eq.  (56)]  of  the  internal  energies,  the  balance 
equation  of  the  gas-mixture  internal  energy 


d^  +  V-(puy)  =  -V-Ju  +  uv 


(59) 
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is  easily  deduced  from  the  summation  of  Eq.  (57)  on  the  subscripts  i,  6.  The  internal- 
energy  diffusive  flux  and  production  on  the  right-hand  side  of  Eq.  (59)  read  respec¬ 
tively 


n  ti 


X!  Juil5 

i=  1  <5=1 

(60) 

n  £i 

yi  yy  Uvjs 

i= 1  <5=1 

(61) 

Equation  (59)  is  not  in  its  final  form.  There  is  more  to  say  about  the  internal-energy 
production  in  consequence  of  the  principle  of  total-energy  conservation.  The  com¬ 
pletion  will  be  done  in  Sec.  7.4. 

Thermal  equilibrium  prevails  when  all  temperatures  equalize  to  a  common 
temperature  T.  This  situation  should  arise  as  a  particular  solution  of  the  multi¬ 
temperature  scheme  embodied  in  Eqs.  (57)  and  (58),  assuming  that  the  component 
internal-energy  diffusive  fluxes  and  productions  are  correctly  prescribed.  An  al¬ 
ternative  approach,  possible  when  there  is  sufficient  (experimental)  evidence  that 
supports  the  idea  as  a  useful  approximation  accurate  enough  to  reflect  realism,  con¬ 
sists  in  the  presumptive  imposition  of  thermal  equilibrium  as  a  shortcut  to  spare  the 
numerical  costs  of  dealing  with  the  mathematical  complexity  of  Eqs.  (57)  and  (58). 
In  this  manner,  the  details  associated  with  Eq.  (57)  are  given  up  and  Eq.  (59)  is  used 
directly  for  the  determination  of  the  gas-mixture  specific  internal  energy  with  the 
provision  that,  now,  a  phenomenological  relation  is  needed  for  the  internal-energy 
diffusive  flux  appearing  on  the  left-hand  side  of  Eq.  (60).  A  phenomenological  re¬ 
lation  for  the  internal-energy  production  appearing  on  the  left-hand  side  of  Eq.  (61) 
is  not  needed  because  its  expression  is  fixed  by  the  imposition  of  total-energy  con¬ 
servation  (see  Sec.  7.4).  Obviously,  the  thermodynamic  relations  (58)  are  still  ap¬ 
plicable  with  Tis  =  T;  thus,  the  temperature  of  the  gas  mixture  follows  from  the 
resolution  of 

n  £i 

PU  =  J2J2  PiUis(T,  Vi)  (62) 

i=  1  <5=1 


7.4  Matter  energy 


According  to  the  considerations  of  Sec.  7.1,  specific  matter  energy  is  defined  as  the 
sum 


=  u  + 


2 


(63) 
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of  kinetic  energy  and  gas-mixture  internal  energy.  The  corresponding  balance  equa¬ 
tion 


+  V  ■  (pemv)  =  -  V  •  (Ju  -  T  ■  v) 

+  uv  —  T  :  Vv  +  pcv  ■  E  -  JQ  •  v  x  B  (64) 

follows,  therefore,  from  the  sum  of  Eqs.  (54)  and  (59).  The  inspection  of  the  right- 
hand  side  of  Eq.  (64)  provides  the  matter-energy  diffusive  flux 

JEm  =  Ju  -  T  ■  V  (65) 


and  production 

em,v  =  uv  —  T  :  Vv  +  pc\  ■  E  -  JQ  ■  v  x  B 


(66) 


In  the  presence  of  the  electromagnetic  field,  one  is  not  entitled  to  assume  the  matter- 
energy  production  as  unconditionally  vanishing.  The  further  addendum  to  account 
for  is  the  production  of  the  electromagnetic  energy 

em,v  =  ~j  ■  E  (67) 


Its  expression  derives  from  arguments  related  to  the  work  done  by  the  electromag¬ 
netic  field  when  electric  charges  are  displaced  within  it  [40].  So,  for  consistency 
with  the  physical  phenomenology,  the  principle  of  total-energy  conservation  must 
be  enforced  as 


+  ~  0 


(68) 


The  substitution  of  Eqs.  (66)  and  (67),  the  latter  expanded  according  to  Eq.  (33), 
into  Eq.  (68)  leads  to  the  following  important,  full  of  physical  significance,  expres¬ 
sion  of  the  internal-energy  production 

uv  =  T  :  Vv  +  JQ  •  (E  +  v  x  B)  (69) 


The  Joule  effect  appears  naturally  in  Eq.  (69)  and  is  represented  by  the  electro¬ 
magnetic  term  linked  exclusively  to  the  conduction-current  density.  With  regard 
to  this  point,  it  seems  worth  mentioning  that  sometimes  the  electromagnetic  term 
on  the  right-hand  side  of  Eq.  (67)  is  erroneously  confused  as  being  responsible 
for  the  Joule  effect.  An  important  conclusion  to  be  drawn  from  Eq.  (69),  with  a 
view  to  Eq.  (61),  is  that  not  only  the  combined  action  of  medium  deformation  and 
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tensional  state  but  also  the  flowing  of  a  conduction  current  through  the  electromag¬ 
netic  field  concurs  to  induce  non-equilibrium  excitation  of  the  molecular  degrees  of 
freedom.  How  the  repartition  of  the  converted  amount  of  energy  takes  place  among 
the  molecular  degrees  of  freedom  can  be  ascertained  only  when  the  expressions  of 
the  productions  uVjis  are  explicitly  known.  In  any  case,  Eqs.  (61)  and  (69)  together 

n 

indicate  that  in  multi-temperature  circumstances,  only  ^  —  1  productions  need 

2=1 

the  assignment  of  phenomenological  relations,  and  that  such  a  necessity  does  not 
exist  in  the  event  of  thermal  equilibrium. 


The  availability  of  Eq.  (69)  leads  to  recast  Eq.  (64)  into  the  final  form 


dpem 

dt 


+  V  •  (pemv) 


— V  •  (Ju  —  X  •  v)  +  j  •  E 


(70) 


and  to  the  completion  of  the  balance  equation  [Eq.  (59)]  of  the  gas-mixture  internal 
energy  which  now  reads 

+  V  ■  H  =  —V  •  Ju  +  T  :  Vv  +  JQ  •  (E  +  v  x  B)  (71) 

at 

Equation  (70)  or  Eq.  (71)  can  replace  anyone  of  Eq.  (57)  in  the  set  of  the  governing 
equations. 


7.5  Electromagnetic  energy 


The  recognition  of  the  existence  of  the  electromagnetic  energy  and  the  derivation  of 
its  associated  balance  equation  are  achieved  by  following  a  procedure  very  similar 
to  the  one  worked  out  for  the  electromagnetic  momentum,  that  is,  through  a  skillful 
transformation  [40,43,46,47]  of  the  electromagnetic-energy  production  [Eq.  (67)] 
by  taking  advantage  of  the  Maxwell  equations  [Eqs.  (36)  and  (37)].  Appendix  D 
provides  the  mathematical  details.  The  final  result  already  cast  in  accordance  with 
Eq.  (14)  reads 

^[|  (E2  +  c2B 2)]  =  -V  •  (e0c2E  xB)-j-E  (72) 

Equation  (72)  indicates  explicitly  that  the  electromagnetic  field  contains  energy  dis¬ 
tributed  in  space  with  density  —  {E2  +  c2B2),  transported  through  space  with  flux 

given  by  the  Poynting  vector  e0c2 E  x  B,  and  exchanged  with  the  energy  of  mat¬ 
ter  with  the  production  rate  —j  •  E.  Once  again,  Eq.  (72)  is  not  an  independent 
equation;  there  is  no  new  physical  information  in  it  that  is  not  already  contained  in 
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the  Maxwell  equations.  The  considerations  made  in  this  regard  with  respect  to  the 
balance  equation  of  electromagnetic  momentum  [Eq.  (52)]  apply  here  unvaried. 

Equation  (72)  reflects  the  structure  of  Eq.  (14)  but  there  have  been  attempts  [13,14] 
to  adapt  it  for  the  purpose  of  fitting  the  structure  of  Eq.  (16).  The  adaptation  is  based 
on  the  transformation  of  electric-field  intensity  and  magnetic  induction  between 
two  reference  systems  in  the  non-relativistic  approximation 


E'  =  E  +  vxB  (73) 

B;  =  B  -  4v  x  E  (74) 

cz 

The  primed  reference  system  is  identified  with  the  one  attached  to  the  generic  fluid 
particle  during  its  motion.  The  basic  step  is  the  evaluation  of  the  Poynting  vector  in 
the  primed  reference  system  [13,14,20] 

£0C2 E'  xB'  =  £0C2 E  x  B  -  j(E2  +  c2B2)\  +  TM  •  v  (75) 

by  taking  advantage  of  the  transformations  (73)  and  (74).  Appendix  E  contains  the 
mathematical  details.  It  is  then  a  simple  matter  to  solve  Eq.  (75)  for  £0c2 E  x  B  and 
to  substitute  the  resulting  expression  into  Eq.  (72)  to  obtain  an  alternative  balance 
equation  of  the  electromagnetic  energy 


|[|(B2  +  c2B2)]  +  V.[|(B2  +  c2B2)v]  = 

-  V  •  (£0C2  E'  xB'-TM-v)-j-E  (76) 

Equation  (76)  reflects  the  structure  of  Eq.  (16)  and  shows  an  interesting  and  re¬ 
markable  similarity  with  Eq.  (70).  From  its  perspective,  electromagnetic  energy 
is  converted  with  matter  and  diffused  with  flux  £„c2E'  x  B'  —  TM  ■  v.  This  view¬ 
point  shares  many  analogies  with  the  one  discussed  in  Sec.  5  in  relation  to  the 
magnetic-induction  equation  [Eq.  (40)].  In  this  case  also,  there  is  a  conceptual  ob¬ 
jection,  already  hinted  at  by  Napolitano  [16],  that  obscures  the  appeal  of  this  adap¬ 
tation  and  of  its  consequent  interpretation.  The  transformations  (73)  and  (74)  are 
rigorously  valid  only  between  two  reference  systems  in  uniform  rectilinear  mo¬ 
tion  with  respect  to  each  other.  Thus,  they  are  not  complete  if  the  primed  refer¬ 
ence  system  is  attached  to  the  generic  fluid  particle  because  the  latter  is  acceler¬ 
ated  (a  =  d\/dt  +  v  •  Vv).  As  explicitly  emphasized  by  Feynman  [40],  transfor¬ 
mations  of  electric-field  intensity  and  magnetic  induction  between  two  reference 
systems  in  relative  accelerated  motion  do  depend  on  the  acceleration.  One  may 
wonder  whether  or  not  the  terms  connected  with  acceleration  that  should  appear  in 
Eqs.  (73)  and  (74)  are  negligible  in  the  non-relativistic  approximation.  Besides  the 
fact  that  general  transformations  including  acceleration  seem  to  be  found  nowhere 
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in  the  literature,  the  question  appears  to  be  a  moot  argument  not  worth  grappling 
with  because  even  if  a  positive  answer  is  found,  one  cannot  reconcile  Eq.  (76)  with 
the  physical  fact  that  the  electromagnetic  field  and  its  properties  momentum  and 
energy  are  attached  to  space. 


7.6  Total  energy 


According  to  Eq.  (68),  all  the  electromagnetic  energy  that  disappears  locally  reap¬ 
pears  as  matter  energy  or  viceversa,  exactly  in  the  same  guise  of  what  happens  to 
momentum.  The  sum  of  Eqs.  (70)  and  (72),  therefore,  provides  the  balance  equa¬ 
tion  of  total  energy 

^[pem  +  | (E2  +  c2 B2)]  =  -  V  •  (pvem  +  J,  -  T  ■  v  +  e0c2E  x  B)  (77) 

For  the  purpose  of  numerical  calculations,  Eq.  (77)  is  perfectly  equivalent  to  either 
Eq.  (70)  or  Eq.  (71)  but,  on  the  contrary  of  the  latter  equations,  it  does  not  present 
any  burdensome  production  term. 


7.7  Mechanisms  of  energy  conversion 


A  summary  of  the  productions  relative  to  kinetic  energy,  internal  energy  and  elec¬ 
tromagnetic  energy  is  illustrated  in  Table  1 .  The  electromagnetic-energy  production 
[Eq.  (67)]  has  been  expanded  according  to  Eq.  (33).  The  tabulation  gives  a  visual 
representation  of  the  possible  mechanisms  of  energy  conversion.  Thus,  electromag¬ 
netic  energy  is  converted  partly  in  kinetic  energy  (pcv  ■  E)  through  the  action  of  the 
electric  field  on  the  convection  current  and  partly  in  internal  energy  (JQ  •  E)  through 
the  action  of  the  electric  field  on  the  conduction  current.  In  turn,  kinetic  energy 
is  converted  in  internal  energy  via  the  interplay  between  medium  deformation  and 
tensional  state  (T  :  Vv),  and  through  the  combined  action  of  the  conduction  cur¬ 
rent  and  the  magnetic  induction  (JQ  ■  v  x  B).  The  Joule  effect  [J0  ■  (E  +  v  x  B)] 
is  the  conjoint  manifestation  of  two  different  conversion  mechanisms  of,  respec¬ 
tively,  electromagnetic  and  kinetic  nature.  A  more  complete  characterization  of  the 
energy-conversion  schematism  illustrated  in  Table  1  covering  aspects  of  reversibil¬ 
ity  and  irreversibility  presupposes  the  explicit  knowledge  of  the  entropy  production. 
By  definition,  the  latter  identifies  the  irreversible  processes  and,  being  subdued 
to  the  second  law  of  thermodynamics  that  guarantees  its  non-negativity,  imposes 
an  inviolable  direction  arrow  on  some  of  the  conversion  pathways  existing  among 
the  corresponding  terms  in  Table  l.The  entropy  production  will  be  dealt  with  in 
Sec.  9.2. 
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Table  1 

Mechanisms  of  energy 

conversion. 

energy  form 

production 

kinetic 

— T  :  Vv 

+  pcv  ■  E 

-JQ  V  x  B 

internal 

+T  :  Vv 

+Jq  ■  E 

+JQ  V  x  B 

electromagnetic 

-pc\  ■  E  -JQ  ■  E 

8  Concluding  considerations  related  to  the  governing  equations  in  open  form 


The  governing  equations  in  open  form  surveyed  in  the  preceeding  sections  embrace 
the  physical  phenomenology  comprising  the  hypersonic  flow  field  of  an  ionized, 
but  not  polarized  and  magnetized,  gas  mixture  and  the  presence  of  the  electromag¬ 
netic  field.  For  quick  reference,  they  are  summarized  in  Tables  2-5  according  to 
several  alternative  but  physically  equivalent  options.  Regardless  of  the  selected  op¬ 
tion,  the  set  of  equations  is  not  operative  yet  because  it  contains  the  variables  requir¬ 
ing  the  assignment  of  phenomenological  relations.  These  variables  identify  the 
fundamental  disciplines  that  converge  into  the  foundational  framework  on  which 
hypersonics  rests,  namely  thermodynamics  chemical  kinetics  (£fc),  thermal 
kinetics  ( uVji$ ),  diffusion  theory  (Jm.,T,  JUi4),  and  call  for  the  selection  of  models 
apt  to  represent  in  an  as  accurate  as  possible  manner  the  physical  behaviour  man¬ 
ifested  by  a  given  real  medium  under  the  specific  circumstances  characteristic  of 
a  given  application.  The  latter  requirement  materializes  through  the  assignment  of 
the  thermodynamic  model  (all  fo)  for  the  gas  mixture  and  of  the  phenomenologi¬ 
cal  relations  establishing  the  link  between  the  unknowns  £k,  uVii§ ,  Jmi,  T,  JUi6  and 
the  basic  unknowns,  and/or  their  gradients,  of  the  flow  field.  Only  then,  the  equa¬ 
tions  in  the  governing  set  acquire  the  prerequisite  closed  form  necessary  to  proceed 
towards  the  achievement  of  their  mathematical  solution.  It  seems  appropriate  at 
this  point  to  emphasize  that  the  seemingly  incomplete  character  of  the  governing 
set  in  open  form  should  not  hinder  at  all  the  development  of  algorithms  for  the 
numerical  solution  of  the  differential  equations  that  belong  to  the  set.  On  the  con¬ 
trary,  such  a  development  is  highly  desirable.  As  a  matter  of  fact,  it  wifi  never  be 
stressed  enough  how  much  convenient  it  is  for  the  efficient  resolution  of  the  flow 
field  that  algorithm-development  studies  would  concentrate  on  the  governing  set  in 
open  form  as  main  target.  This  is  a  proposition  that  certainly  implies  an  ambitious 
programme  but  the  prospective  benefits  are  too  appealing  to  be  ignored  and  the 
idea  to  be  hurriedly  dismissed.  If  such  a  programme  succeeds  then  the  phenom¬ 
enological  relations  become  relegated  to  the  role  of  subroutines ,  interchangeable 
according  to  the  specific  necessities  of  a  given  application,  and  the  architecture  of 
the  numerical  kernel  will  feature  the  extraordinary  useful  flexibility  of  being  gen¬ 
erally  applicable  and  independent  from  the  specific  physical  behaviour  of  the  real 
medium. 
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Table  2 


Governing  equations  relative  to  gas-mixture  composition 


equations 

Eq. 

unknowns 

no.  of  equations 
or  sub.  range 

dp  r  ■ 

+  V  •  (pi\)  =  —V  ■  Jin,  +  ^  £,kvkiMi 

k= 1 

(17) 

Pi 

i  =  1,  •  ■  • ,  n  —  s 

f)  * 

-g  +  V-  (p*v)  =  -V  .  J„. 

(22) 

* 

Pi 

3  =  !,•••,  s 

n  1 

Jm*  =  55  TF^miai0Mj 

i= 1  M' 

(23) 

Jm* 

3  =  !,•••,  s 

0?i 

II 

<-i.  * 

EMI 

£1$ 

q 

iS>. 

(27) 

Pi 

3  = 

pi  =  E  p*i  &ji 

3= 1 

(29) 

Pi 

i  =  n  —  s+1,  •  •  • ,  n 

|  +  V.(pv)=0 

(21) 

P 

1 

n 

Pc  =  ~ eiVA  Y  —&is 

Mi 

(13) 

Pc 

1 

n  l 

Jq  =  —  6j^a55  ~j\/T^mi<7is 

i= 1  1 

(31) 

Jq 

3 

3  =  pc\  +  Jq 

(33) 

3 

3 

The  construction  of  the  thermodynamic  model  is  a  task  belonging  to  the  realm  of 
statistical  thermodynamics  (see  Sec.  7.3).  There  are  a  few  options  available  con¬ 
cerning  the  derivation  of  the  phenomenological  relations  for  the  unknown  produc¬ 
tions  and  diffusive  fluxes.  One  can  seek  recourse  to  irreversible  thermodynamics 
[20-24,35-37,67-70],  to  the  more  sophisticated  kinetic  theory  of  gases  [71-77]  or 
to  experimental  investigation.  In  practice,  the  phenomenological  relations  emerge 
as  the  outcome  of  a  concerted  effort  involving  all  three  options  to  different  degrees 
of  depth.  The  approach  relying  on  irreversible  thermodynamics  is  preferable  to  get 
started  in  the  derivation  endeavour  because,  although  its  findings  may  have  some¬ 
times  narrow  limits  of  validity  from  a  quantitative  point  of  view,  it  proceeds  in  a 
conceptually  straightforward  manner  from  the  exploitation  of  the  entropy  produc¬ 
tion  and  of  the  second  law  of  thermodynamics,  it  is  not  affected  (not  to  the  same 
extent,  at  least)  by  the  overwhelming  mathematical  cumbersomeness  and  complex- 
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Table  3 


Governing  equations  relative  to  momentum 


equations 

Eq. 

unknowns 

no.  of  equations 

+  V  ■  (pv  v)  =  V  •  T  +  pcE  +  j  x  B 

(48) 

V 

3 

or 

-^(pv  +  e0E  x  B)  =  -V  •  (pvv  -  T  -  XM) 

(53) 

V 

3 

Tm  =  £„(EE  -  l-E2 U)  +  £0c2(BB  -  l-B2 U) 

(51) 

Tm 

6 

ity  of  the  detailed  kinetic  theories  and,  above  all,  it  offers  a  depth  of  insight  that 
goes  a  long  way  in  the  direction  of  understanding  the  transport  processes  at  work 
in  the  flow  field  and  of  recognizing  the  associated  driving  forces. 


9  Linear  irreversible  thermodynamics 

9.1  Preliminary  remarks 


The  linear  theory  of  irreversible  thermodynamics  will  be  revisited  in  the  follow¬ 
ing  sections  in  conformity  with  the  prescription  of  thermal  equilibrium.  The  as¬ 
sumption  that  thermal  equilibrium  prevails  among  the  molecular  degrees  of  free¬ 
dom  of  the  components  is  a  recurrent  characteristic  shared  by  authors  that  follow 
the  irreversible-thermodynamics  approach.  Some  [23,37]  even  go  further  and  as¬ 
sume  mechanical  equilibrium.  Thermal  equilibrium  implies  the  possibility  of  deal¬ 
ing  with  one  single  temperature  and,  obviously,  brings  in  great  simplification;  on 
the  other  hand,  it  restricts  the  applicability  domain  of  the  ensuing  phenomenolog¬ 
ical  relations.  The  sole  attempts  the  present  author  is  aware  of  that  ventured  into 
a  thermal-disequilibrium  analysis  were  made  by  Woods  [22],  Napolitano  [35,36], 
and  Morro  and  Romeo  [78-80].  However,  the  treatments  proposed  by  Woods  and 
by  Napolitano  share  similarities  that  contain  elements,  bearing  on  the  definition  of 
the  driving  forces  connected  with  the  occurrence  of  multiple  temperatures,  appar¬ 
ently  not  yet  completely  freed  from  conceptual  objections.  Similarly,  Morro  and 
Romeo  did  not  consider  the  internal  structure  of  the  molecules;  in  other  words, 
they  implicitly  assumed  for  each  component  the  thermal  equilibrium  among  its 
molecular  degrees  of  freedom.  More  work  is  certainly  needed  to  improve  knowl¬ 
edge  in  this  department  of  irreversible  thermodynamics.  The  motivation  justifying 
the  choice  adopted  here  resides  mainly  in  the  intention  to  put  the  emphasis  on 
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Table  4 


Governing  equations  relative  to  energy 


equations 

Eq. 

unknowns 

no.  of  equations 
or  sub.  range 

+  V  •  (PiuiSv)  =  -  V  •  JUi5  +  uv,i6 

(57) 

'U'iS 

<5=1, ■ 
i  =  1,  • 

. .  f. 

i 

•  • ,  n 

uiS  —  T- iS  ^  dTs  )  ~  uiS\TiS,Vi) 

(58) 

Tis 

<5=1, ■ 
i  =  1,  ■ 

■■Ji 
•  • ,  n 

Vi  =  1/  Pi 

Vi 

i  =  1,  ■ 

•  • ,  n 

n  £i 

pu  =  T,J2  PiUis 

i= 1  <5=1 

(56) 

u 

1 

V2 

em  =  u  +  — 

(63) 

1 

or 


—  V  •  JL,  +  T  :  Vv  +  JQ  •  (E  +  v  x  B) 

(71) 

u 

1 

n  li 

Ju  =  Jui5 

(60) 

Ju 

3 

i=l  <5=1 

d  QiVbifo 

Qt  +  ^  ■  ( PiUiSV )  —  V  •  Ju-j  +  UV}i$ 

(57) 

^i5 

n 

£«<- 

2—1 

1 

n  li 

PU  =  y]y]  PiuiS 

(56) 

'U'iS 

1 

i= 1  <5=1 

rr2  fdfis/Tis\  _ [rr  ^ 

7); ^  ^  <9T,5  /  ui&\TiSi  vi) 

(58) 

<5  =  1, 

*  =  !>••■ 

■,4 

,n 

Vi  =  I/Pi 

Vi 

*  =  !,••■ 

,n 

V2 

em  =  «  +  y 

(63) 

1 

the  peculiarities  of  the  physical  phenomenology  connected  with  the  existence  of 
the  electromagnetic  field  and,  for  that  purpose,  to  keep  the  mathematical  analysis 
relieved  from  tangential  or,  even,  unnecessary  complexity.  Nevertheless,  the  multi¬ 
temperature  phenomenology  is  quantitatively  important  in  hypersonic  applications 
and  should  not  be  forgotten. 
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Table  4 
Continued 


equations 

Eq. 

unknowns 

no.  of  equations 
or  sub.  range 

or 

dpem  ,  . 

at  +  V'(',e"v,= 

-  V  •  (Ju  -  X  •  v)  +j  E 

(70) 

1 

n  U 

Ju  =  ^  Juii 

(60) 

Ju 

3 

i= 1  <5=1 

+  V  •  (fHUisv)  =  -  V  •  JLU  +  uv,iS 

(57) 

'U'iS 

n 

E«(-i 

i=  1 

n  £i 

Pu  =  y]y]  Piui8 

(56) 

'U'iS 

l 

i=  1  <5=1 

_  rr2  fdfis/TiS\  _ {rr  ^ 

uiS  ^  QT g  J  ui5\Ti8ivi) 

(58) 

Tis 

S  =  l,---,£i 
i  =  1,  •  •  • ,  n 

Vi  =  1/ Pi 

Vi 

i  =  1, ■ ■ ■ ,  n 

V2 

em  =  u  +  — 

(63) 

u 

1 

or 

^  [Pe  m  +  y  (-®2  +  c2-82)]  = 

-V  •  (/9vem  +  Ju  -  X  •  v  +  e0c2E  x  B) 

(77) 

1 

n  li 

Ju  =  ^  ^  y  ^  Ju,;a 

(60) 

Ju 

3 

i=  1  5=1 

+  V  •  (piUisV)  =  -V  •  Ju,4  + 

(57) 

'U'iS 

n 

Ee (-1 

i=l 

n  t?; 

pu=y]y]  Piui& 

(56) 

'U'iS 

l 

i= 1  5=1 

rr2  f  ^fis/Tis\  ___  (rr  ^ 

(58) 

TiS 

S  =  l,---,£i 
i  =  1, ■ ■  • ,  n 

Vi  =  1/ Pi 

Vi 

i  =  1, • ■  ■ ,  n 

V2 

em  =  It  +  y 

(63) 

u 

1 
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Table  5 

Governing  equations  relative  to  the  electromagnetic  field 


equations 

Eq. 

unknowns 

no.  of  equations 

dB  „  „ 

—  +  V  x  E  =  0 
ot 

(36) 

B 

3 

<9E  ™ 

Eo^rr  -  EoC  V  X  B  =  -j 

Ot 

(37) 

E 

3 

or 

1  d2(j)  2  pc 

c2  8E  V^+£0 

(44) 

<t> 

1 

1  d2  A  2  1  . 

c2  dt2  “  V  A  +  £oC2  7 

(45) 

A 

3 

B  =  V  x  A 

(43) 

B 

3 

-4*  l+j 

1 

-e- 

> 

I 

II 

(42) 

E 

3 

9.2  Entropy-balance  equation  and  entropy  production 


The  thermodynamics  of  the  irreversible  processes  hinges  on  the  entropy  balance 
equation 

^+V-(pv4  =  -V-  Js  +  4  (78) 

Equation  (78)  represents  a  formal  balance  involving  entropy  density  4,  diffusive 
flux  Js,  and  production  4„,  and  acquires  operational  character  only  after  that  entropy 
diffusive  flux  and  production  are  given  more  explicit  expressions.  This  is  achieved 
by  taking  advantage  of  the  thermodynamic  fundamental  relation  [24,39,68,81]  of 
the  gas  mixture  in  the  entropic  formulation  which,  under  the  assumption  of  thermal 
equilibrium,  reads  formally 

J=  J(u,v,aua2:. . .  ,an)  (79) 


In  Eq.  (79),  v  =  1/p  is  the  specific  volum  of  the  gas  mixture.  In  view  of  the  forth¬ 
coming  analysis,  it  is  appropriate  to  recall  that  the  first  derivatives  of  the  function 
(79)  represent  the  state  equations  in  the  entropic  scheme  and  define,  respectively, 
temperature,  pressure  and  chemical  potentials 
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1 

T 

P_ 

T 

T 


v,a 


s,a 


i  =  1,  • 


’  5  n 


u,v,aj:£i 


(80) 

(81) 

(82) 


of  the  gas  mixture.  Also,  for  reasons  of  simplicity  in  the  notation,  it  is  convenient  to 
work  out  the  mathematical  details  with  the  aid  of  the  substantial-derivative  operator 
D  / Dt  =  d/dt  +  \-V,  whose  intervention  in  the  balance  equations  stems  from  the 
equivalence  [recall  also  Eq.  (16)] 


dpg 

dt 


+  V • (p\g) 


(83) 


Equation  (83)  is  a  direct  consequence  of  total-mass  conservation  [Eq.  (21)].  Taking 
into  account  the  definitions  (80)-(82),  the  application  of  the  substantial-derivative 
operator  to  the  function  (79)  leads  to 


Dj  1  Du 

P15i  =  TP~Dt 


p  Dv  JC  Pi  Dat 

TP~Dt  ~^TP~Dt 
2=1 


(84) 


The  steps  to  achieve  the  explicit  form  of  the  entropy-balance  equation  appear  now 
evident.  The  first  step  consists  in  obtaining  the  volume-balance  equation.  The  quick¬ 
est  way  is  to  set  g  =  v  in  Eq.  (83)  and  to  take  into  account  that  pv  —  1;  hence 


Dv  _ 


(85) 


Counterintuitively  as  it  may  appear  at  first  sight,  the  term  on  the  right-hand  side  of 
Eq.  (85)  does  not  represent  the  divergence  of  the  diffusive  flux  of  volume  because, 
by  its  very  definition,  a  diffusive  flux  cannot  depend  on  the  velocity  vector  [see 
Eq.  (15)];  this  inescapable  requirement  suggests  that  a  volume -production  (yv  = 
V  ■  v)  interpretation  is  the  sole  possible  for  the  term  in  question.  An  alternative 
explanation  reaching  the  same  conclusion  can  be  found  in  Ref.  [19].  The  second 
step  involves  the  substitution  of  the  component-mass  balance  equations  [Eq.  (17)], 
the  internal-energy  balance  equation  [Eq.  (71)]  and  the  volume-balance  equation 
[Eq.  (85)]  into  Eq.  (84),  and  the  consequent  rearrangement  to  reflect  the  structure 
of  Eq.  (78).  This  brings  to  the  identification  of  the  entropy  diffusive  flux 

Js  =  =  Ju  -  E  ^  (86) 

1  i= i  1 


and  of  the  entropy  production 
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4  =  :  Vv 

1  n  i  n  i 

-  ^(Ju  -  E^JE  •  VT  -  -  53  Jmi  ■  V/0  +  ^Jq  •  (E  +  V  X  B) 

1  i=  1  J  i=l  1 

1  ^ 

+  |V'v-^E^A.  (87) 

k=  1 

The  affinities  of  the  chemical  reactions  are  defined  as  linear  combinations 

n 

Ak  =  E  "kiMifii  (88) 

i=  1 

of  the  chemical  potentials.  The  entropy  production  [Eq.  (87)]  is  not  in  its  final  form 
but  must  be  furtherly  refined  for  the  purpose  of  the  present  context.  The  refinement 
consists  in  performing  the  standard  expansion  of  the  stress  and  velocity-gradient 
tensors  to  bring  forth  their  isotropic,  traceless-symmetric  (or  deviatoric)  and  anti¬ 
symmetric  parts.  This  operation  presupposes  some  familiarity  with  tensor  analysis; 
extensive  help,  in  this  regard,  is  provided  in  Refs.  [22,23,37,71,74].  The  stress  ten¬ 
sor  separates  in  the  sum 

T  =  TtU  +  Tos  (89) 


of  the  isotropic  part,  which  contains  the  (scalar)  mean  normal  stress  built  as  the 
arithmetic  average  (7t  =  1/3  T  :  U)  of  the  diagonal  terms  of  the  stress  tensor,  and  of 
the  traceless-symmetric  part  T\  The  antisymmetric  part  of  T  vanishes  identically 
because  the  stress  tensor  remains  symmetric  in  the  absence  of  polarization  and 
magnetization.  The  same  expansion  performed  on  the  velocity  gradient  yields 

Vv=  ^(V-v)U  +  (Vv):  +  (Vv)a  (90) 

The  isotropic  part  [1/3(V  ■  v)U]  describes  volume  deformation  of  fluid  particles, 
the  traceless  symmetric  part  [(Vv)0s]  describes  their  shape  deformation,  and  the  an¬ 
tisymmetric  part  [(Vv) a]  describes  their  rigid  rotation  with  1/2  V  x  v  as  angular 
velocity.  With  the  help  of  relatively  simple  tensor  algebra,  it  can  be  shown  that,  tak¬ 
ing  into  account  Eqs.  (89)  and  (90),  the  first  term  on  the  right-hand  side  of  Eq.  (87) 
becomes 

T  :  Vv  =  7t(V  ■  v)  +  Tos  :  (Vv)os  (91) 


The  substitution  of  Eq.  (91)  and  of  the  expression  of  the  conduction-current  density 
in  terms  of  the  component-mass  diffusive  fluxes  [Eq.  (31)]  into  Eq.  (87)  leads  to 
the  desired  final  form  of  the  entropy  production 
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4  =  :  (Vv): 

1  n  -in 

-  ~(J„  -  ■  VT  -  -  5:  J„,  .  [Vft  +  eA^E  +  v  x  B)] 

z=l  i= 1  1 

-^[-(tt  +  p)]V-v-  ~J2£kAk  (92) 

k= 1 

In  irreversible-thermodynamics  parlance,  Eq.  (92)  shows  the  generalized  fluxes  and 
the  generalized  forces  that  produce  entropy  and,  in  so  doing,  identifies  the  processes 
responsible  for  promoting  irreversibility  in  the  flow  field.  In  addition  to  the  habit¬ 
ual  effects,  namely  shape  [(Vv)Ds]  and  volume  (V  •  v)  deformation  of  the  medium, 
nonequilibrium  chemical  reactions  (Ak  0),  thermal  (VT)  and  chemical-potential 
(V Hi)  disuniformities,  also  the  electromagnetic  field  enlists  in  the  category  of  the 
irreversible  processes  through  contributions  that  add  to  the  chemical-potential  gra¬ 
dients  and  connect  to  the  mass  diffusive  fluxes  of  the  electrically  charged  compo¬ 
nents.  Generalized  fluxes  and  forces  are  summarized  in  Table  6. 


Table  6 

Generalized  fluxes  and  forces  contributing  to  the  entropy  production  [Eq.  (92)] 


flux 

force 

tensorial  order 

subscript  range 

(Vy): 

2 

1  n 
m(Ju 

i=  1 

VT 

1 

Jm; 

Vq,  +  eNAj^j-( E  +  v  x  B) 

1 

i  =  1,  •  •  • ,  n 

-(n  +  p) 

V  -v 

0 

6c 

Ak 

0 

k  =  1,  •  •  ■ ,  r 

The  core  of  irreversible  thermodynamics  is  the  exploitation  of  the  expression  of  the 
entropy  production  [Eq.  (92)]  and  of  the  second  law  of  thermodynamics 

4  >  0  (93) 


which  endorses  its  unconditional  non-negativity  as  an  ineluctable  fact  of  nature,  for 
the  purpose  of  obtaining  suitable  expressions  of  the  generalized  fluxes  in  terms  of 
the  generalized  forces.  In  this  regard,  the  validity  of  the  Curie  postulate  is  assumed 
here.  The  postulate  affirms  that  generalized  fluxes  of  a  given  tensorial  order  can 
depend  only  on  generalized  forces  of  the  same  tensorial  order.  This  is  a  simplifica¬ 
tion  because  the  postulate  is  rigorously  applicable  only  for  isotropic  media;  that  is 


1  -40 


RTO-EN-AVT-162 


Hypersonic-Flow  Governing  Equations  with  Electromagnetic  Fields 


not  the  case  in  the  presence  of  the  electromagnetic  field  because  the  existence  of 
the  magnetic  induction  removes  the  medium  isotropy.  Thus,  the  assumption  of  the 
postulate  validity  implies  the  negligibility  of  the  couplings  among  each  generalized 
flux  and  all  the  generalized  forces  of  different  tensorial  order.  The  motivation  be¬ 
hind  this  choice  is  the  same  adduced  at  the  end  of  Sec.  9.1  to  justify  the  assumption 
of  thermal  equilibrium,  namely  the  intention  to  emphasize  the  essential  characteris¬ 
tics  due  to  the  presence  of  the  electromagnetic  field  and,  at  the  same  time,  minimize 
mathematical  complexity.  Readers  interested  in  the  general  case  that  includes  the 
neglected  cross  effects  are  referred  to  the  specialized  literature  mentioned  in  Sec.  8 
(Refs.  [22,23]  in  particular).  To  compensate  to  some  extent  for  the  limitations  en¬ 
suing  from  the  Curie  postulate,  the  anisotropy  imposed  on  the  medium  by  the  mag¬ 
netic  induction  will  be  taken  into  account  in  the  tensorial  nature  of  the  transport 
coefficients  [22,23,7 1 ,74,75,77,82-90] . 

In  compliance  with  the  invoked  Curie  postulate,  it  is  convenient  to  emphasize  the 
subdivision  of  the  entropy-production  terms  with  same  tensorial  order  by  rephras¬ 
ing  Eq.  (92)  as 

4;  4,2  f  -4; _ |  T  dvo  (94) 


in  which 


4  —  — Ts  • 

*Ai,2  ji  Vo  ' 


Vv) 


(95) 


4,1  —  rp2  (Ju  AfiJm,)  ■  VT 
1  i= 1 


1  n 

yEJ»-[Vft+elVA-|(E  +  vxB)] 

i=l  1 


1  1  r  • 

•f, o  =  —  ^  [— (^  T  p)]  V  ■  v  —  —  £,k-Ak 

'  k=i 


(96) 

(97) 


Given  the  independence  of  the  three  groups  (95),  (96)  and  (97),  the  non-negativity 
property  of  the  entropy  production  [Eq.  (93)]  is  passed  on  to  each  of  them.  They 
will  be  dealt  with  separately  in  the  following  sections. 


9.3  Tensorial  (second-order)  generalized  force 


The  traceless  symmetric  part  of  the  velocity  gradient  is  the  sole  generalized  force 
of  second  tensorial  order  that,  conjointly  with  the  traceless  symmetric  part  of  the 
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stress  tensor,  participates  in  the  entropy  production  [Eq.  (95)] 


(98) 


The  application  of  the  Curie  postulate  leads  to  the  relatively  simple  phenomeno¬ 
logical  relation  [22,23,74] 


(99) 


that  resembles  the  more  familiar  Newton  law.  The  dynamic  viscosity  u.  in  Eq.  (99) 
is  a  fourth-order  tensor  due  to  the  anisotropy  introduced  by  the  magnetic  induction. 
Its  determination  is  beyond  the  reach  of  irreversible  thermodynamics  and  it  must  be 
considered  an  import  either  from  the  kinetic  theory  of  gases  or  from  experiments. 
A  thorough  analysis  of  the  dynamic-viscosity  tensor  is  beyond  the  purpose  of  the 
present  context.  Here,  it  is  sufficient  to  recall  that  most  of  the  81  scalar  coefficients 
required  to  compose  the  tensor  turn  out  to  vanish  and  that  its  complete  characteri¬ 
zation  is  limited  to  the  knowledge  of  only  five  scalar  coefficients.  The  sign  of  some 
of  them  must  be  necessarily  positive  for  compliance  with  the  non-negativity  of  4, 2. 
References  [22,23,74]  provided  detailed  descriptions  of  structure  and  properties  of 
the  dynamic-viscosity  tensor. 

9.4  Vectorial  generalized  forces 

Phenomenological  relations  for  vectorial  generalized  forces  arise  from  the  entropy 
production  related  to  the  first  tensorial  order 


1 


dv,\  —  rp2  (Ju  /-fiJm,- )  ■  VT 


2—1 


(100) 


In  this  regard,  one  has  to  confront  a  multifaceted  situation  because  the  combination 
of  generalized  fluxes  and  forces  appearing  in  Eq.  (100)  is  not  unique.  Different, 
although  equivalent,  forms  of  4,i  exist  in  consequence  of  transformations  to  which 
fluxes  and  forces  can  be  subjected,  a  peculiarity  due  to  the  occurrence  of  mass  dif¬ 
fusion  and  recognized  long  time  ago.  The  equivalence  of  the  several  forms  stems 
from  the  expectable,  and  provable,  invariance  of  the  phenomenological  relations 
for  internal-energy  diffusive  flux  and  component-mass  diffusive  fluxes  with  respect 
to  the  mentioned  transformations.  On  the  other  hand,  the  multiplicity  of  forms  im¬ 
plies  the  non-uniqueness  of  the  heat-flux  definition.  This  fact  should  not  come  as 
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a  surprise  or  be  perceived  as  the  manifestation  of  some  inconsistency  in  the  theory 
because,  as  appropriately  remarked  by  Napolitano  [35,36]  in  1969,  the 

. . .  arbitrariness  in  the  definition  of  heat  flux  stems  from  the  fact  that  any  time 
there  is  a  flux  of  matter  . . .  there  is  also  a  flux  of  energy  . . .  and  it  is  only  a  mat¬ 
ter  of  convention  (or  of  convenience)  to  single  out  one  or  more  of  these  fluxes  of 
energy  and  call  them  heat  flux.  . . .  Any  possible  definition  obviously  leaves  all 
physical  results  unchanged  and  one  should  only  pay  attention  in  using  it  consis¬ 
tently,  i.e.  in  associating  to  it  the  pertinent  expression  for  the  entropy  production 


For  a  given  selection  of  generalized  forces,  the  heat  flux  is  identified  with  the  gen¬ 
eralized  flux  associated  with  the  temperature  gradient.  Thus 

n 

Jq  =  Ju  -  /A-Jm,  =  TJS  (101) 

1=1 

is  the  heat-flux  definition  compatible  with  the  generalized  forces  that  appear  in 
Eq.  (100).  Equation  (101)  is  somewhat  reminiscent  of  the  classical-thermodynamics 
relation  dQ  =  TdS  but  it  has  received  scant  attention  and  its  implications  on  the 
derivation  of  the  phenomenological  relations  have  rarely  been  pursued  in  detail 
[23,24,68].  Another  form  of  4,;1  encountered  in  the  literature  [21,23,35-37]  is  ob¬ 
tained  by  contracting  the  terms  with  the  component-mass  diffusive  fluxes 

•  vr  -  7j„,,  •  d02) 

in  Eq.  (100).  In  this  way,  the  latter  becomes 

1  1  n  1 1  st 

4,i  =  -^Ju  ■  VT  -  -  £  Jm,  •  [TV^  +  eNA^-( E  +  v  x  B)]  (103) 

2—1  1 

and  yields  the  heat-flux  definition 

n 

Jq  =  Ju  =  TJS  +  ^  /i?;  J  m,  (104) 

i= 1 


The  overwhelmingly  preferred  form  [23,35-37]  of  4m,  however,  is  the  one  arrived 
at  from  the  assumption  according  to  which  the  components  of  the  gas  mixture  be¬ 
have  as  perfect  gases.  Indeed,  in  this  case  the  chemical  potential  can  be  expressed 
in  terms  of  specific  enthalpy  h%  and  entropy  4,;  of  the  7-th  component 

fii  =  hi-TJi  (105) 
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and  its  gradient  depends  on  those  of  temperature  and  partial  pressure  pi 


V/ii  =  -J.VT  +  ViVpi 


(106) 


The  substitution  of  Eqs.  (105)  and  (106)  into  Eq.  (100)  leads  to  the  form 


4>,l  —  (Ju  E  hjJmi)  •  VT 


1 


i= 1 


1  n 

-  f  E  J...  '  [»iVPi  +  eN^  (E  +  v  x  B)] 

i= 1  1 


(107) 


with  the  corresponding  heat- flux  definition 


n 


n 


J"'  =  Ju-E  ^ Jnu  =  TJS  -  T  E  4Jm, 


(108) 


The  instances  considered  so  far  by  no  means  exhaust  the  list  of  possibilities.  Other 
choices  [37]  of  the  generalized  fluxes  and  forces,  together  with  the  identification  of 
the  correspondig  heat  flux,  are  possible.  On  the  other  hand,  it  is  neither  necessary 
nor  intended  in  the  present  context  to  embark  in  an  exhaustive  survey.  The  three 
described  forms  of  4M  suffice  to  portray  the  manifoldness  of  the  heat-flux  concept. 
The  irrelevance  of  its  arbitrariness,  pointed  out  in  Napolitano’s  remark,  becomes 
mathematically  evident  and  comprehensible  by  overturning  Eq.  (101),  Eq.  (104) 
and  Eq.  (108)  as 


n 


n 


Ju = j; + e  m.,, = j; = j.7 + e  w« 


(109) 


and  by  keeping  in  mind  that  the  internal-energy  diffusive  flux  is  invariant  with 
respect  to  any  possible  selection  of  generalized  fluxes,  forces  and  corresponding 
heat-flux  definition.  Indeed,  it  is  the  internal-energy  diffusive  flux  to  be  the  physical 
flux  that  really  matters  in  order  to  evaluate  the  thermal  energy  transferred  from  the 
flow  field  into  a  body  or  viceversa. 

The  derivation  of  the  phenomenological  relations  of  first  tensorial  order  pursued 
here  differs  from  the  habitual  treatments  found  in  the  literature.  It  proceeds  along  a 
pathway  that,  starting  from  a  slight  adaptation  of  Eq.  (107),  establishes  a  strong  par¬ 
allelism  with  the  Chapman-Enskog  method  [71-77,82,84-88]  of  the  gas-kinetics 
theory  and  enormously  facilitates  the  import  of  the  transport  coefficients  provided 
by  that  method.  The  transformation  that  takes  care  of  bringing  4>,i  into  the  appropri¬ 
ate  form  presupposes  the  temporary  replacement  of  the  component-mass  diffusive 


1  -44 


RTO-EN-AVT-162 


Hypersonic-Flow  Governing  Equations  with  Electromagnetic  Fields 


flux  in  the  term 


Jmi  •  [viVpi  +  eNAjjL(  E  +  vxB)] 

(HO) 

with  the  intermediate  appearance  of  the  diffusion  velocity  defined 

in  Eq.  (18).  After 

substituting  the  latter  into  Eq.  (110),  taking  into  account  that  f)tv, 
for  brevity 

=  1,  and  setting 

=  Vp*  +  eNAp~(E  +  v  x  B) 

(HD 

Eq.  (107)  can  be  rephrased  as 

n  n 

T4,tl  =  -(Ju  -  E  i)  •  V  InT  -  E  Wi  •  Xi 

2=1  2=1 

(112) 

The  advantage  of  using  Eq.  (112)  rather  than  Eq.  (107)  is  that  the  former  leads  to 
operate  with  the  central  ingredients  [71,72,74]  of  the  Chapman-Enskog  method, 
namely  the  vectors  w*,  X,  and  VlnT.  The  application  of  the  Curie  postulate  to 

Eq.  (112)  provides  the  phenomenological  relations 

n 

W i  =  -J2  '  Xj  -  LiT  ■  V  In  T 

3= 1 

(113) 

n  n 

Ju  -  E  hi  Jms  =  -  E  L«  •  x,  -  LTT  ■  V  In  T 

2=1  2=1 

(114) 

in  which  the  phenomenological  coefficients  L,:j ,  L)T.  LX(  and  Lxx 

are  second-order 

tensors  due  to  the  presence  of  the  magnetic  induction.  Equations  (113)  and  (114) 
mirror  exactly  the  expressions  intervening  in  the  Chapman-Enskog  method  and,  in 
consequence  of  that,  the  tensors  L?J- ,  L,x.  LX)  and  Lxx  can  be  directly  and  explicitly 

related  as 

.  1  pMiMjn 

13  pPi  TP 

(115) 

L'  =  —  DT 

Pi  1 

(116) 

LTj  =  ib'(-B) 

Pi 

(117) 

Lxx  —  Tyl 

(118) 
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to  the  tensors  of  diffusion  DtJ,  of  thermal  diffusion  D(  and  to  the  primed  thermal- 
conductivity  tensor  >,' .  The  prime  attached  to  the  thermal-conductivity  tensor  in 
Eq.  (118)  is  important:  it  represents  the  definition  of  thermal  conductivity  consis¬ 
tent  with  the  scalar  one  that  does  not  include  the  thermodiffusion  ratios  [72,74]. 
The  latter  tensors  are  styled  according  to  the  (scalar)  notation  of  Hirschfelder,  Cur¬ 
tiss  and  Bird  [72]  only  for  reasons  of  convenience  related  to  flow  calculations. 
Yet,  the  diffusion  tensors  Di?  in  Eq.  (115)  must  be  consistent  with  the  conven¬ 
tion  adopted  by  Chapman  and  Cowling  (Eqs.  (18.2,9)  and  (18.3,4)  in  Ref.  [71]} 
and  Ferziger  and  Kaper  (Eq.  (6.3-24)  in  Ref.  [74]}  to  fix  the  indeterminacy  of  the 
diffusion-related  coefficients  appearing  in  the  first-order  perturbation  of  the  distri¬ 
bution  function  in  order  to  make  Eq.  (113)  operational.  The  diffusion  tensors  con¬ 
sistent  with  the  indeterminacy-fixing  convention  adopted  by  Hirschfelder,  Curtiss 
and  Bird  (Eq.  (7.3-31)  in  Ref.  [72]}  are  incompatible  with  Eq.  (113)  because  of  the 
presence  in  it  of  the  vectors  Xj  rather  than  the  vectors  (Eq.  (18.2,6)  in  Ref.  [71]} 


Vxj  +  (Xj  -  % V  lnp  -  pkFk) 


(119) 


P  PP  “i 


habitually  employed  in  kinetic  theory  and  based  mainly  on  the  gradients  of  the 
molar  fractions  xr  In  the  present  case,  the  external  forces  in  Eq.  (119)  arise  from 
the  electromagnetic  field 


(120) 


Another  important  reason  that  justifies  this  choice  was  expressed  by  Ferziger  and 
Kaper  [74] 

We  remark  that  there  is  a  considerable  variation  among  authors  in  the  nomen¬ 
clature  and  definition  of  the  multicomponent  diffusion  and  thermal  diffusion  co¬ 
efficients.  The  definition  [adopted  by  Ferziger  and  Kaper,  as  well  as  Chapman 
and  Cowling]  of  the  multicomponent  diffusion  coefficients  is  consistent  with 
Onsager’s  reciprocity  relations  of  irreversible  thermodynamics  ...  This  consis¬ 
tency  is  of  particular  importance  if  one  desires  to  generalize  the  kinetic  theory 
of  gas  mixtures  to  systems  in  which  the  number  density  of  each  species  is  not 
a  conserved  quantity  —  e.g.,  in  gas  mixtures  in  which  chemical  reactions  occur 
between  constituents  or  in  polyatomic  gases  in  which  transition  occur  between 
states  with  different  internal  degrees  of  freedom  ...  It  has  been  shown  ...  that,  in 
these  cases,  it  is  of  importance  to  choose  the  definition  of  the  transport  coeffi¬ 
cients  in  such  a  way  that  they  are  consistent  with  Onsager’s  reciprocity  relations 
and  that  the  definition  of  the  multicomponent  diffusion  coefficients  given  [here] 
is  therefore  to  be  preferred  to,  e.g.,  the  definition  given  by  Hirschfelder,  Curtiss 
and  Bird  ...  Recently,  J.  Curtiss  ...  has  supported  this  point  of  view. 
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In  Eq.  (115),  M  is  the  average  molar  mass 


M 


-l 


EM*Xi 

i= 1 


(121) 


of  the  gas  mixture.  Equation  (117)  is  an  Onsager  reciprocal  relation.  The  tilde  in¬ 
dicates  matrix  transposition  and  the  attached  functional  dependence  reminds  the 
necessary  inversion  of  the  sign  of  the  magnetic  induction.  Similar  reciprocal  re¬ 
lations  apply  also  to  the  diffusion  and  thermal-conductivity  tensors.  However,  al¬ 
though  important,  reciprocity  details  are  not  of  specific  relevance  in  the  present 
context;  interested  readers  are  referred  to  the  specialized  literature,  particularly  the 
textbooks  by  Woods  [22]  and  De  Groot  and  Mazur  [23]. 


The  phenomenological  relations 


Jmi  pM^M 

t-  3=\ 


.  T 


T 


DJ  •  VT 


(122) 


for  the  component-mass  diffusive  fluxes  follow  from  the  substitution  of  Eqs.  (115) 
and  (116)  into  Eq.  (113)  and  of  the  resulting  diffusion  velocities  into  Eq.  (18).  The 
kinetic-theory  result 


n 

EM*Do=°  (123) 

i= 1 

n 

EDI=°  (124) 

i= 1 

makes  the  condition  (19)  identically  satisfied.  In  the  same  guise,  the  substitution 
of  Eqs.  (117)  and  (118)  into  Eq.  (1 14)  yields  the  phenomenological  relation  for  the 
internal-energy  diffusive  flux 

n  n  I 

Ju  =  E  Mm,  *  E  -Dl  •  Xi  -  X'  •  VT  (125) 

i=  1  i=l  Pi 

In  Eq.  (125),  the  functional  dependence  (— B)  has  been  dropped  from  the  trans¬ 
posed  thermal-diffusion  tensor  in  order  to  simplify  the  notation;  nevertheless,  this 
peculiarity  should  always  be  kept  in  mind. 

In  order  to  derive  phenomenological  relations  analogous  to  Eqs.  (122)  and  (125), 
several  authors  [23,35-37]  followed  an  alternative  procedure  starting  directly  from 
Eq.  (107).  In  their  analyses,  the  term  v^pi  appears  as  the  isothermal  gradient  of 
the  chemical  potential  whose  definition  descends  from  the  following  rearrangement 

ViVpi  =  V/i,  +  JtVT  =  Vj/m  (126) 
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of  Eq.  (106).  Then,  the  generalized  forces  associated  with  the  component-mass 
diffusive  fluxes  in  Eq.  (107)  appear  as 


J^V^  +  eJV^E  +  vxB) 


(127) 


The  ensuing  phenomenological  relations  for  the  component-mass  diffusive  fluxes 
and  the  internal-energy  diffusive  flux  are  further  manipulated  to  bring  forth  the  gra¬ 
dients  of  mass  fractions  (Fick  law)  and  pressure.  The  drawback  of  this  procedure 
resides  in  the  fact  that  the  intervening  tensorial  phenomenological  coefficients,  say 
l_E,  L't,  I_7,  L(t,  are  not  expressible  in  terms  of  the  kinetic-theory  tensors  D  ‘ ,  >,' 
as  simply  as  in  Eqs.  ( 1 15)— (1 1 8)  but,  in  order  to  obtain  the  required  correlation, 
necessitate  an  elaborate  manipulation  involving  the  resolution  of  cumbersome  al¬ 
gebraic  systems  of  equations. 

The  phenomenological  relations  (122)  and  (125)  are  operational  equations  that,  to¬ 
gether  with  Eq.  (Ill)  and  with  adequate  input  from  the  kinetic  theory,  qualify  to 
complement  the  set  of  governing  equations  in  open  form  for  the  purpose  of  nu¬ 
merical  calculations.  However,  noticeable  insight  towards  the  understanding  of  the 
transport  processes  can  be  gained  if  the  role  played  by  the  contributions  compos¬ 
ing  the  generalized  force  X,;  is  made  explicit.  The  substitution  of  Eq.  (Ill)  into 
Eq.  (122)  yields 


with  the  tensorial  coefficients 


(129) 


Equation  (128)  shows  explicitly  how  the  electromagnetic  field  participates  to  the 
build-up  of  the  component-mass  diffusive  fluxes  at  the  same  footing  of  the  partial- 
pressure  and  temperature  gradients.  The  appearance  of  the  mass-fraction  gradients 
is  one  step  away  from  Eq.  (128)  and  is  produced  by  expanding  the  partial  pressures 
with  the  aid  of  the  perfect-gas  relations  p^  =  Xip  and  xt  =  arMt/M .  Accordingly, 
Eq.  (128)  becomes 


k= 1 


with  the  tensorial  coefficients 
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D/fc  = 


Mk 


Mn 


n 


XE^Evt?)  P 

3= 1 


M 


(Djfc  — 


Din) 


pm 

’  i 


1  M  n 
, ;  y  Dj  f); 
P  M  j=i 


(131) 

(132) 


Equation  (130)  warns  against  the  uncareful  practice  of  limiting  the  expressions  of 
the  component-mass  diffusive  fluxes  to  the  sole  Fick  law,  namely  the  first  term  on 
the  right-hand  side.  Even  in  the  event  of  negligible  thermodiffusion  and  vanishing 
electromagnetic  field,  that  widely  assumed  law  may  not  be  sufficient  to  charac¬ 
terize  completely  the  component-mass  diffusive  fluxes  because  the  tensorial  coef¬ 
ficients  P“,  being  constructed  according  to  combinations  of  the  diffusion  tensors 
very  similar  to  those  that  appear  in  Eq.  (131),  could  have  magnitudes  compara¬ 
ble  to  the  tensorial  coefficients  T>ik.  Thus,  it  seems  safer  not  to  discard  a  priori 
the  pressure-gradient  contribution  in  Eq.  (130)  but  to  let  its  negligibility,  if  any, 
arise  spontaneously  during  the  course  of  numerical  simulations.  Equation  (131)  il¬ 
lustrates  the  sometimes  forgotten  fact  that  the  diffusion  coefficients  appearing  in 
Fick  law  must  be  obtained  from  appropriate  adaptation  of  those  provided  by  the 
Chapman-Enskog  method. 


The  substitution  of  Eq.  (128)  into  Eq.  (31)  leads  to  a  more  explicit  and  meaningful 
expression  of  the  conduction-current  density 

n 

Jq  =  E  Kj  •  Vp?  +  Xe  ■  (E  +  v  x  B)  +  Xj  •  VT  (133) 

3  = 1 


and  to  the  recognition  of  the  pressoelectrical-c onductivity  tensors 

*  i= 1 


(134) 


the  electrical-conductivity  tensor 


Xc  —  (eiVA)  pj^/pi  E  E 


(135) 


in  rigorous  agreement  with  the  expression  derived  by  Chapman  and  Cowling  [71], 
and  the  thermoelectrical-conductivity  tensor 


-I  71 

XT  =  eivEy  ^  DJ 


T  Mi 
1=1  1 


(136) 


Equation  (133)  evidences  unambiguously  that  the  generalized  Ohm  law  alone  may 
not  be  sufficient  to  characterize  the  conduction-current  density  in  an  ionized  gas 
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mixture  because,  even  in  the  event  of  negligible  thermodiffusion,  the  contributions 
associated  with  the  partial-pressure  gradients  may  not  be  negligible  with  respect  to 
the  diffusion  promoted  by  the  electromagnetic  field.  These  considerations  have  a 
negative  impact  on  the  applicability  of  the  magnetic-induction  equation  discussed 
in  Sec.  5.  At  the  same  time,  Eq.  (135)  casts  skepticism  about  the  physical  mean¬ 
ingfulness  of  (actually,  it  severely  admonishes  against)  undertaking  flow-field  sim¬ 
ulations  in  the  presence  of  the  electromagnetic  field  in  which  the  gas  mixture  is 
assimilated  to  a  perfect  gas  without  electrically  charged  carriers  but  still  with  a 
finite  (scalar)  electrical  conductivity.  This  is  a  flagrant  infringement  of  physical 
coherence  because  the  electrically-conductivity  tensor  vanishes  identically  in  the 
absence  of  electrically  charged  components,  as  Eq.  (135)  unconditionally  indicates 

( ®is  °js  —  0)- 

An  expression  structured  in  the  same  style  of  Eq.  (128)  can  be  obtained  also  for  the 
internal-energy  diffusive  flux  by  substituting  Eqs.  (Ill)  and  (128)  into  Eq.  (125); 
the  final  outcome  reads 

n 

JU  =  E  Pi  •  +  Fu  ■  (E  +  v  x  B)  -  X"  •  VT  (137) 

3= 1 

with  the  tensorial  coefficients 


p]=£(%-D,  +  ^ 


i=  1 


liT  p ,  MjMj 

M2 


Pi 


P 


Do) 


(138) 


P  =  eiVA^ 


i=l 


_  T  (Tis  p,  M, 


P  '  M2 


E  D  ijPj 


a 


(139) 


1  n 

x;  =  X'  +  -E^°I  (140) 

1  i= 1 

There  are  two  aspects  related  to  Eq.  (137),  usually  not  encountered  in  traditional 
hypersonic-flow  theory,  that  may  have  implications  for  spacecraft  thermal  protec¬ 
tion  during  planetary  reentry.  As  clearly  shown  in  Eq.  (137),  there  is  a  direct  con¬ 
tribution  arising  just  from  the  mere  existence  of  the  electromagnetic  field.  In  par¬ 
ticular,  the  influence  of  the  magnetic  induction  disappears  on  the  wall  (v  =  0)  of  a 
reentry  body  but  the  influence  of  the  electric  field  could  still  play  a  role.  Moreover, 
Eq.  (137)  indicates  also  that,  on  the  contrary  of  the  widespread  habitude  acquired 
within  the  realm  of  a  scalar  transport  theory,  the  component  of  J,  normal  to  the 
wall  of  a  reentry  body  depends  not  only  on  the  normal  gradients  of  partial  pres¬ 
sures  and  temperature  but  also  on  their  tangential  gradients  [23].  With  the  goal  in 
mind  of  achieving  a  controlled  mitigation  of  the  internal-energy  diffusive  flux,  this 
occurrence  gives  cause  for  the  necessity  to  produce  an  accurate  assessment  of  the 
importance  of  the  off-diagonal  terms  in  the  transport  tensors.  More  specifically, 
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a  thorough  understanding  of  the  anisotropic  character  of  the  transport  tensors  un¬ 
der  the  influence  of  the  electromagnetic  field  appears  to  be  a  worthwhile  target. 
Chapman  and  Cowling  [71]  pioneered  the  study  of  this  subject  but  the  cases  they 
considered  were  subdued  to  simplifications  (steady  states,  slight  or  full  ionization, 
binary  mixtures,  absence  of  chemical  reactions,  mixtures  macroscopically  at  rest, 
and  the  like)  which  appear  (sometimes  very)  restrictive  when  looked  at  from  the 
perspective  of  hypersonics  applications.  Chapman  and  Cowling’s  findings  and  con¬ 
clusions  have  reverberated  throughout  the  literature  [11,22,74];  unfortunately,  the 
accompanying  simplifications  seem  to  have  been  somehow  diluted  during  the  rever¬ 
beration  process.  Thus,  central  issues  are  still  in  demand  of  attention.  For  example, 
it  is  important  to  identify  the  circumstances  according  to  which  the  anisotropy  pos¬ 
sessed  by  the  transport  tensors  is  or  is  not  negligible,  to  recognize  if  the  anisotropy 
works  against  or  in  favor  of  the  sought  mitigation,  to  investigate  the  existence  of 
potential  ways  to  enhance  the  anisotropy  effects  and  drive  them  towards  the  sought 
mitigation.  The  work  already  done  [71,74,75,77,82-90]  requires  improvement  in 
the  generality  of  the  findings  and,  if  necessary,  in  their  adaptation  to  routine  use  in 
numerical  calculations. 


9.5  Scalar  generalized  forces 


The  derivation  of  the  phenomenological  relations  relative  to  the  scalar  generalized 
forces  is  based  on  the  entropy  production  related  to  the  zero  tensorial  order 

4,o  =  -^[-(tc  +  p)]V-v  -  ^J2ikAk  >  0  (141) 

k=  1 


The  application  of  the  Curie  postulate  to  Eq.  (141)  provides  the  phenomenological 
relations 


r 


^  ^  k^k  ZvvV  *  V 

k= 1 

r 

(142) 

^  v  lkj*Aj  I'kv^^  *  V 

3= 1 

(143) 

The  scalar  phenomenological  coefficients  lyk,  lyy,  lkj,  lky  form  an  [(r  +  1)  x  (r  + 
1)]  matrix  which  must  be  positive  definite  in  order  to  have  the  inequality  (141) 
identically  satisfied.  Onsager  reciprocal  relations  exist  also  for  these  coefficients; 
plenty  of  information  is  available  in  the  specialized  literature.  The  scalar  phenom¬ 
enological  coefficients  appearing  in  Eqs.  (142)  and  (143)  have  not  been  investi¬ 
gated  with  the  same  attention  devoted  to  the  kinetic -theory  tensors  encountered 
in  Secs.  9.3  and  9.4  because  it  was  soon  recognized  [23,35,36,69]  that  the  linear 
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combinations  of  the  affinities  in  Eqs.  (142)  and  (143)  have  very  narrow  limits  of 
validity.  The  coefficient  Zvv  is  better  known  as  the  bulk- viscosity  coefficient 

ZVV  =  u„  (144) 


and  is  the  only  one  which  has  received  attention  because  it  characterizes  the  direct 
effect  of  the  medium  compressibility  (V  ■  v)  on  the  generation  of  normal  stresses. 
Even  with  the  limitation  imposed  by  their  restricted  applicability,  Eqs.  (142)  and 
(143)  retain  importance  from  a  theoretical  point  of  view  because  they  reveal  unam¬ 
biguously  the  existence  of  cross  effects  between  chemical  reactions  and  medium 
compressibility  [21,35-37].  This  aspect  has  been  systematically  ignored  by  au¬ 
thors  concerned  with  CFD  applications.  Notwithstanding  a  few  daring  attempts 
[91-93]  towards  theoretical  progress  on  the  subject  matter,  the  concern  expressed 
by  Napolitano  [35,36]  in  1969  that 

Not  much  is  known  nor  is  done  ...  to  ascertain  the  order  of  magnitude  of  this 
type  of  coupling  . . . 

is,  unfortunately,  still  applicable  today. 

The  disregard  of  cross  effects  brings  Eq.  (142)  into  the  simplified  form 


7l  -  —p  —  u,,V  ■  v 


(145) 


which  finds  widespread  use  in  the  applications,  most  often  accompanied  by  the 
ulterior  approximation  of  vanishing  bulk  viscosity  (u„  =  0).  In  the  same  guise, 
Eq.  (143)  would  become 

ik  =  —hk-Ak  (146) 


However,  chemical-reactions  rates  are  never  used  in  the  form  (146)  but,  rather,  in 
the  non-linear  form 
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(147) 


that  features  explicit  dependence  on  the  component  partial  densities  and  that  some¬ 
times  is  referred  to  as  the  law  of  mass  action  [23,35,36].  Equation  (147)  arises  from 
a  blend  of  theoretical  considerations  and  (mainly)  experimental  evidence  [76,94- 
99].  The  reaction  constants  K/k.  K,\  are  in  reality  functions  of  the  temperature,  cus¬ 
tomarily  assumed  in  accordance  with  a  modified  version 

ftf'b  =  CTve~Ea/KsT  (148) 
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of  the  Arrhenius  law,  and  are  related  to  each  other  through  the  chemical-equilibrium 
constant  defined  in  terms  of  the  component  concentrations 


(149) 


A  theoretical  justification  for  Eq.  (148)  can  be  produced  only  in  particularly  sim¬ 
ple  cases.  In  it,  KB  is  the  Boltzmann  constant  (1.3806503TCT23  JK_1)  while  the 
constant  C,  the  temperature  exponent  77,  and  the  activation  energy  Ea  are  derived, 
in  practice,  by  utilizing  that  expression  as  an  interpolation  formula  constructed  on 
experimentally  determined  data  and,  then,  collected  in  databases  [100,101]  readily 
available  to  serve  the  needs  of  the  applications.  The  reliability  of  Eq.  (148)  becomes 
weak  in  the  high-temperature  range  due  to  serious  difficulties  encountered  in  carry¬ 
ing  out  accurate  measurements.  The  uncertainties  of  the  reaction  constants  at  high 
temperatures  is  a  well  known  and  long  standing  problem  in  demand  of  resolution. 
In  the  chemical-equilibrium  limit  [(Ah  — ►  0)  =>  Ak/RaTc  1],  Eq.  (147)  reduces 
[35-37,69]  to  Eq.  (146)  with  the  phenomenological  coefficient  given  by 


(150) 


The  subscript  “e”  attached  to  the  partial  densities  in  Eq.  (150)  indicates  chemical- 
equilibrium  values  corresponding  to  a  local  couple  of  thermodynamic  parameters 
[(T,  v)  or  (T,p)];  R0  is  the  universal  gas  constant  (8.314472  J-K”1). 

Notwithstanding  the  widespread  use  of  Eq.  (147)  in  CFD,  it  is  important  to  keep 
in  mind  that  the  premises  on  which  it  is  constructed  reflect  the  chemical  kinetics 
of  a  perfect-gas  mixture  contained  in  a  constant  volume,  macroscopically  at  rest, 
without  spatial  disuniformities,  and  whose  population  distributions  over  the  quan¬ 
tum  energy  states  associated  with  the  molecular  degrees  of  freedom  can  still  be 
sufficiently  assimilated  to  a  Boltzmann  distribution  [35,36].  The  chemical-kinetics 
phenomenology  described  by  Eq.  (147)  is,  therefore,  substantially  simpler  than  the 
more  complex  one  accompanying  hypersonics  situations.  For  example,  some  un¬ 
resolved,  but  of  extreme  relevance,  issues  are  the  failure  of  Eq.  (147)  to  bring  into 
account  the  influence  of  the  medium  compressibility,  the  intrinsic  difficulty  to  adapt 
it  to  multitemperature  situations,  and  its  incompatibility  with  quantum-state  popu¬ 
lations  that  do  not  follow  a  Boltzmann  distribution. 


10  Conclusions 

The  understanding  of  the  coupling  between  the  hypersonic  flow  field  of  an  ionized 
gas  mixture  and  the  electromagnetic  field  is  conceptually  straightforward  in  the  ab- 
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sence  of  polarization  and  magnetization  effects.  The  recognition  of  the  weaknesses 
affecting  the  assumptions  on  which  the  derivation  of  the  largely  used  magnetic- 
induction  equation  is  based  indicates  that  the  adoption  of  either  directly  the  Maxwell 
equations  [Eqs.  (36)  and  (37)]  or  the  field  equations  for  the  electromagnetic  scalar 
and  vectorial  potentials  [Eqs.  (44)  and  (45)]  represents  a  better  choice  reflecting  a 
more  consistent  physical  description  of  the  electromagnetic  field.  The  cost  to  pay  is 
the  development  of  appropriate  numerical  algorithms  that,  especially  in  the  case  in 
which  the  field  equations  of  the  potentials  are  chosen,  may  turn  out  to  be  substan¬ 
tially  different  from  those  in  use  since  many  years  in  computational  fluid  dynam¬ 
ics.  Obviously,  the  coverage  of  the  physical  phenomenology  relative  to  polarization 
and  magnetization  is  desirable  for  reasons  of  both  scientific  interest  and  engineer¬ 
ing  necessity.  Thus,  the  open-form  governing  equations  should  be  extended  for  that 
purpose.  The  extension  should  be  accomplished  in  a  bilateral  manner,  that  is,  not 
only  by  using  the  appropriate  Maxwell  equations  but  also  by  introducing  the  re¬ 
quired  ingredients  on  the  fluid-dynamics  side  such  as,  for  example,  the  specific 
angular  momentum,  its  associated  energy,  their  respective  balance  equations,  and 
so  on.  Subsequently,  the  formalization  of  state-to-state  thermal  kinetics  and  radia¬ 
tion  together  with  their  integration  in  the  open-form  governing  equations  represent 
a  wishful  endeavour  towards  completion. 

The  open-form  governing  equations  must  be  complemented  with  phenomenolog¬ 
ical  relations.  In  this  respect,  the  method  of  irreversible  thermodynamics,  based 
on  the  analysis  of  the  entropy  production,  proves  useful  to  recognize  generalized 
fluxes  and  generalized  forces  at  work  in  the  flow,  achieves  deep  insight  into  the 
physics  of  their  cross  couplings,  and  provides  linear  phenomenological  relations 
that,  although  with  the  exception  of  the  chemical-reaction  rates,  go  rather  far  in 
satisfying  the  needs  of  engineering  applications.  The  method  straightforwardly  re¬ 
veals  how  the  influence  exerted  by  the  electromagnetic  field  on  the  diffusion  of 
(electrically-charged)  component  masses  and  internal  energy  becomes  manifest. 
There  is  a  novelty  with  respect  to  more  habitual  fluid-dynamics  situations  without 
electromagnetic  field:  the  presence  of  the  magnetic  induction  destroys  the  isotropy 
of  the  flow.  The  introduced  anisotropy  calls  for  particular  attention  to  the  possibility 
of  cross  coupling  between  generalized  fluxes  and  forces  of  different  tensorial  order 
and  to  understand  the  importance  of  the  tensorial  nature  of  the  transport  coefficients 
with  the  goal  in  mind  of  exercising  control  over  the  distributions  of  internal-energy 
diffusive  flux  imposed  on  the  surfaces  of  a  body  during  planetary  reentry.  Multi¬ 
temperature  circumstances  should  not  be  underestimated.  There  have  been  efforts 
to  bring  the  multitemperature  phenomenology  within  the  reach  of  irreversible  ther¬ 
modynamics  but  the  actual  status  of  the  theory  demands  further  improvements. 
This  appears  a  mandatory  task  to  pursue  in  order  to  incorporate  correctly  the  ef¬ 
fects  originating  from  the  thermal  disequilibrium  among  the  molecular  degrees  of 
freedom  of  the  components. 

The  final  goal  of  the  theoretical  fluid  dynamicist  is  the  construction  of  an  ade¬ 
quate  mathematical  apparatus  with  the  proper  physics  built  in  it.  According  to  the 
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previous  considerations,  there  is  plenty  of  work  in  the  theoretical  arena  awaiting 
investigation  and  understanding.  New  problems  join  old  ones,  some  of  which  were 
identified  and  formulated  long  time  ago  but  are  still  unresolved  today.  In  this  re¬ 
spect,  it  can  be  safely  maintained  that  two  antithetical  approaches  emerge.  The 
minimalist  approach  consists  in  throwing  into  the  equations  a  simplified  physics, 
most  of  the  times  even  stretched  beyond  its  domain  of  applicability,  and  rushing 
to  the  computer  in  the  hope  that  the  discarded  physics  is  unimportant  and  that  ei¬ 
ther  numerically  predicted  flows  match  experimental  evidence  within  an  acceptable 
approximation  or  sensitivity  analyses  produce  conservative  results.  This  is  a  risky 
approach.  For  one  thing,  accurate  experimental  evidence  in  hypersonic  regime  may 
be  difficult  to  gather;  in  addition,  the  necessity  imposed  by  engineering  applications 
to  produce  results  within  awfully  short  time  limits  induces  to  yield  at  the  temptation 
of  numerically  fine-tuning  the  physical  models,  an  exercise  that  sometimes  is  cam¬ 
ouflaged  under  the  epithet  of  validation.  The  maximalist  approach  is  characterized 
by  an  attitude  of  thoroughness  towards  the  physics  translated  into  the  equations; 
conceptual  difficulties  are  not  avoided  but  are  confronted  without  hesitation  in  or¬ 
der  to  push  further  ahead  the  comprehension  of  the  phenomenology  dealt  with. 
This  is  a  lengthy  approach,  certainly  incompatible  with  the  shortness  of  the  time 
limits  mentioned  above.  On  the  other  hand,  it  is  inescapable  when  the  minimalist 
approach  leads  into  a  situation  of  stagnating  failure.  The  choice  between  the  two 
approaches  is  a  matter  of  personal  taste.  The  last  years  have  seen  an  unquestionable 
unbalance  towards  the  former  approach  sponsored  by  an  unjustified  rush  to  compu¬ 
tational  fluid  dynamics.  Probably,  the  most  convenient  way  to  go  about  is  the  one 
that  takes  advantage  of  both.  When  dealing  with  new  physical  complexities,  the 
minimalist  approach  is  acceptable  to  start  with  but  resources  should  also  be  allo¬ 
cated  to  promote  a  parallel  line  of  development  following  the  maximalist-approach 
philosophy.  In  the  author’s  opinion,  the  latter  is  a  sine  qua  non  condition  to  achieve 
progress  in  such  a  complex  discipline  as  hypersonics  is. 
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Appendices 


A  Magnetic-induction  equation 

Neglecting  the  displacement-current  density  e0dE/dt,  Eq.  (37)  becomes 

J=V2VxB  (A.l) 

Expanding  the  electric-current  density  in  Eq.  (A.l)  according  to  Eq.  (33)  gives  the 
conduction-current  density 

JQ  =  £0c2 V  xB-pcv  (A. 2) 

The  presumed  applicability  of  the  generalized  Ohm  law  [Eq.  (38)]  brings  Eq.  (A. 2) 
into  the  form 

Xc  (E  +  v  x  B)  =  £0c2V  x  B  —  pcv  (A. 3) 

from  which  the  electric  field 

2 

E  =  xB--v-vxB  (A. 4) 

X,  xe 

can  be  obtained.  The  substitution  of  Eq.  (A. 4)  into  Eq.  (36)  and  the  subsequent 
expansion  of  the  terms  affected  by  the  curl  operator  leads  to  Eq.  (40). 


B  Field  equations  for  scalar  and  vector  potentials 


Taking  into  account  that  the  divergence  of  a  curl  and  the  curl  of  a  gradient  give 
always  a  null  result,  Eq.  (35)  implies  that  the  magnetic  induction  can  be  expressed 
in  terms  of  a  vector  potential  and  of  an  arbitrary  function  of  space  and  time  as 

B  =  Vx(A'  +  Vf)  (B.l) 


The  replacement  of  three  scalar  quantities,  namely  the  components  of  the  vector  B, 
with  four  scalar  quantities,  namely  the  components  of  the  vector  A'  and  the  scalar 
function  0,  introduces  an  arbitrariness  about  the  vector  potential  A  =  A'  +  V0 
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which,  in  order  to  be  resolved,  requires  the  assignment  of  a  conveniently  chosen 
scalar  condition.  With  this  provision,  Eq.  (B.l)  can  be  simply  rephrased  as 


B  =  V  x  A 


(B.2) 


The  field  equation  for  the  vector  potential  derives  from  Eq.  (37).  Performing  the 
curl  on  both  sides  of  Eq.  (B.2)  yields 


VxB  =  Vx(VxA)  =  V(V  ■  A)  -  V2A 


(B.3) 


The  substitution  of  Eq.  (B.2)  into  Eq.  (36)  provides  the  form 

V  x  E  +  -^-V  x  A  =  0 

at 


(B.4) 


which,  after  permutation  of  the  time  derivative  with  the  gradient  operator,  turns  into 

<9  A. 

V  x  (E  +  — )  =  0  (B.5) 


The  quantity  in  brackets  is,  thus,  irrotational  and  can  be  expressed  as  the  gradient 
of  a  scalar  potential 


dA 

E+dt=  -v< 


(B.6) 


From  Eq.  (B.6),  the  electric-field  intensity  is  obtained  in  terms  of  the  potentials  as 

E  =  -v^-nr  (B.7) 

at 


The  substitution  of  Eqs.  (B.3)  and  (B.7)  into  Eq.  (37)  leads  to 


£0C2[V(V  ■  A)  -  V2A]  =j-  e0  + 


r)A , 

dt  ' 


(B.8) 


which  can  be  rearranged  as 

1  d2A 


1  dq ) 


c2  dt2 


-V2A  +  V(V-A+--^)  = 


c2  dt  £0C2" 


(B.9) 


Inspection  of  Eq.  (B.9)  suggests  that  the  gauge  transformation 

1  <90 

+VA  =  ° 


(B.10) 
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rather  than  the  gauge  V  ■  A  =  0,  customarily  imposed  in  magnetostatics,  is  the 
most  convenient  choice  to  resolve  the  arbitrariness  about  A.  The  condition  (B.10)  is 
known  as  the  Lorentz  gauge  [40,47].  Following  it,  Eq.  (B.9)  reduces  to  the  simpler 
Eq.  (45). 

The  field  equation  for  the  scalar  potential  derives  from  Eq.  (34).  The  substitution 
of  Eq.  (B.7)  into  Eq.  (34)  gives 

-V  •  (V0)  -  V  •  ^  =  —  (B .11) 

at  e0 

Equation  (B  .11)  is  equivalent  to 

-vv-4v'a= -  <b-12) 

at  en 


The  term  V  ■  A  is  obtained  from  the  Lorentz  gauge  [Eq.  (B.10)] 


V  A 


1  <90 
c2  dt 


and  substituted  into  Eq.  (B.12)  to  obtain  Eq.  (44). 


(B.13) 


C  Transformation  of  the  electromagnetic  body  force  and  derivation  of  the 
balance  equation  of  electromagnetic  momentum 


The  electric-charge  density  and  the  electric-current  density  can  be  obtained  from 
Eqs.  (34)  and  (37)  as 


pc  =  £0V  •  E 


(C.l) 


3 


o„  <9E 

£0C  V  X  B  -  £0  — 
at 


(C.2) 


and  substituted  into  Eq.  (49)  to  yield 
gv  =  £0EV  ■  E  +  e0c2(V  x  B)  x  B 


<9E 

~dt 


£0—  X  B 


(C.3) 


At  the  same  time,  one  can  take  advantage  of  Eqs.  (35)  and  (36)  to  construct  the 
identity 


0 


£0c2BV  •  B  +  £0(V  x  E  + 


x  E 


(C.4) 
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and  to  expand  it  into  the  expression 


<9B 


0  =  £0c  BV  •  B  +  £0(V  x  E)  x  E  +  £0—  x  E 


(C.5) 


somewhat  symmetrical  with  respect  to  Eq.  (C.3).  The  sum  of  Eqs.  (C.3)  and  (C.5) 
yields 


kv  =  £0c2[BV  •  B  +  (V  x  B)  x  B] 


d 


+  £o[EV  •  E  +  (V  x  E)  x  E]  -  —  (£0E  x  B) 

Taking  into  account  that  the  terms  in  squared  brackets  can  be  contracted  as 


(C.6) 


V  •  (BB-  -B2U) 
v  2  ’ 

(C.7) 

V  •  (EE--E2U) 
v  2  ' 

(C.8) 

Eq.  (C.6)  can  be  rephrased  in  the  final  form 

&,  =  V  •  [£0c2(BB  -  l-B2 U)  +  £o(EE  -  l-E2 U)]  -  ^(£0E  x  B) 


(C.9) 


which,  considering  the  definitions  (49)  and  (51),  coincides  with  Eq.  (50). 


D  Transformation  of  the  electromagnetic-energy  production  and  derivation 
of  the  balance  equation  of  electromagnetic  energy 


Equation  (37)  can  be  resolved  to  express  the  electric-current  density  as 

2__  „  <9E 

J  =  £0C  V  X  B  -  So¬ 
ot 


(D.l) 


Taking  into  account  Eq.  (D.l),  the  electromagnetic-energy  production  becomes 

9  <9E  9  dE2/  2 

3  ■  E  =  £0c2E  -V  x  B  —  £„E  •  —  =  £0c“E  V  x  B  -  £0 — (D.2) 

at  at 


The  first  term  on  the  right-hand  side  of  Eq.  (D.2)  can  be  easily  transformed 

E  •  V  x  B  =  —V  •  (E  x  B)  +  B  •  V  x  E  (D.3) 
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to  bring  forth  the  curl  of  the  electric-field  intensity.  The  substitution  of  Eq.  (36)  into 
Eq.  (D.3)  yields 

nt»  a  jd'2  /o 

E  •  V  x  B  =  -V  •  (E  x  B)  -  B  •  —  =  -V  •  (E  x  B) - -1—  (D.4) 

K  J  dt  K  ’  dt 


so  that  Eq.  (D.2)  finally  takes  the  form 

j  E  =  -V  (e0c2E  x  B)  -  ^[| (E2  +  c2B 2)]  (D.5) 

that  coincides  with  Eq.  (72). 


E  Poynting- vector  transformation 

From  Eqs.  (73)  and  (74),  the  vector  product  between  E'  and  B'  can  be  expanded  as 
follows 

E'xB'  =  ExB  — 1E  xfvxE)-Bx(vxB) 
c2 

-^(vxB)  x  (vxE)  (E.l) 

c2 

The  second  term  on  the  right-hand  side  of  Eq.  (E.  1)  expands  as 

E  x  (v  x  E)  =  (E  ■  E)v  -  (E  ■  v)E  =  {E2 U  -  EE)  ■  v 

E2  E2 

=  — v+(  — U-EE)-v  (E.2) 

Similarly,  the  third  term  on  the  right-hand  side  of  Eq.  (E.l)  expands  into 

r> 2  td2 

B  x  (v  x  B)  =  —v  +  (— U  -  BB)  v  (E.3) 

The  fourth  term  on  the  right-hand  side  of  Eq.  (E.l)  is  a  little  bit  more  complicated. 
It  expands  as 

(v  x  B)  x  (v  x  E)  =  (v  x  B)  •  Ev  =  E  •  (v  x  B)v  =  v  •  (B  x  E)v 

=  -v  •  (E  x  B)v  (E.4) 
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In  the  non-relativistic  approximation,  therefore,  the  term 

^r(v  x  B)  x  (v  x  E)  =  — \\  •  (E  x  B)v  (E.5) 

c2  cz 

is  negligible  with  respect  to  E  x  B  and  can  be  disregarded  in  Eq.  (E.l).  Taking  into 
account  the  definition  [Eq.  (51)]  of  the  Maxwell  stress  tensor,  the  substitution  of 
Eqs.  (E.2)  and  (E.3)  into  Eq.  (E.l)  and  the  subsequent  multiplication  by  the  factor 
£0c2  leads  to  Eq.  (75). 
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